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Introduction and Background 

Monitoring and understanding crustal movements has been identified as a 
major goal of NASA's Office of Space and Terrestrial Applications, with im- 
plications for earthquake prediction, seismic risk evaluation, engineering 
application, resource exploration, and crustal dynamics in general. In de- 
veloping this new space based capability it is essential to all aspects of 
the program (independent of which system is being considered) to incorporate 
and build upon currently available information regarding crustal dynamics. 

In the earliest stages of development, the successful deployment of any 
crustal movement monitoring network utilizing space techniques will depend 
upon obtaining accurate information on such questions as: where should land 
based stations be deployed to most effectively monitor crustal movements 
in both stable areas for calibration purposes and active areas where move- 
ments are most likely to be detected? What precision is necessary to de- 
tect crustal movements in different areas? How often should movements be 
repeated (this is particularly important in light of the possible complex 
temporal behavior of crustal movements)? In addition, when attempting to 
measure continent wide deformation and plate tectonic movements, regional 
and local movements (whether due to tectonic or non-tectonic causes) will 
be a source of noise which could greatly complicate interpretation of the 
observations. Tha ability to identify and remove such effects could prove 
to be the fundamental limit on the applicability of space techniques to 
studies of plate tectonic movements and continent wide deformation. 

For the past 2 years a group at Cornell University under NASA sponsor- 
ship has been using U.S. releveling data to investigate these questions. 

This report briefly summarizes the most recent results of this effort. 
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Recent Results (for details of particular studies see appropriate appendix) 

I. Neotectonic Deformation, Near Surface Movements and Systematic Errors 

in U.S. Releveling Measurements; Implications for Earthquake Prediction 

(Appendix II, III) 

Analyses of U.S. releveling mearsurements indicate that derivative crustal 
movement estimates may reflect tectonic deformation, near-surface movements, 
and/or systematic errors. Discriminating the contributions of these factors 
is especially crucial for unambiguous geodetic detection of possible precursory 
seismic deformations. IVhile reliable leveling measurements of co-seismic and 
post-seismic movements are well documented for some of the larger (M > 6) dip- 
slip earthquakes, leveling evidence for pre-seismic motion is generally sparse 
and often ambiguous. Subtle earthquake-related motions may be masked by both 
aseismic movements and systematic errors. For example, deep magma injection 
and surficial groundwater withdrawal are two mechanisms which have been docu- 
mented to cause surface movements which, under some circumstances, could be 
misidentified as seismic-related. Of more concern, perhaps, are systematic 
measurement errors. Topography dependent errors are an exceptionally trouble- 
some type, perhaps affecting as much as 20% of U.S. leveling. However, other 
varieties of systematic error also contribute to the uncertainty. Discrepan- 
cies between leveling and tide gauge data and within nets of leveling alone 
suggest large, long baseline accumulations of error. In many cases, aseismic 
and erroneous contributions cannot be unequivocally determined ex post facto . 
However, a comprehensive examination of the NGS crustal movement data base, 
representing a large sampling of the entire U.S. Level Net, provides perspec- 
tive and criteria needed to begin to recognize movement directly related to 
earthquake activity. 


Perhaps the most extensive set of relevant measurements exists in 
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southern California, where much attention has recently been focused. 

Reevaluation of some of these leveling observations indicates that while 
some appear to reflect tectonic deformation, others are suspect because of 
indications of systematic errors and/or near-surface, non-tectonic movements. 
Specifically, possible preseismic movements reported for the 1971 San 
Fernando earthquake in the vicinity of the earthquake fault as well as approxi- 
mately 30 km northwest of the epicenter may be due to systematic errors. 
Movements near the San Gabriel fault, initially ascribed to the Palmdale Bulge 
and more recently to preseismic effects of the 1971 San Fernando earthquake 
apparently reflect near-surface sediment compaction due to water table fluc- 
tuations. Similarly, there is strong evidence of contamination by rod cali- 
bration errors in the releveling observations used to define the southwestern 
portion of the "Palmdale Bulge" (Llano to Azusa, California). The reality 
of the "Palmdale Bulge" itself must be questioned in view of this reevaluation. 

In contrast, possible tilting southwest of Palmdale between 1961 and 1964 is 
not easily related to systematic errors or near-surface movements and thus may 
represent tectonic deformation. Whether this tilt anomaly was due to preseismic 
effects of the San Fernando earthquake or mechanically separate tectonic 
event is presently unknown. 

II . Vertical Crustal Movements in the Vicinity of the 1951 Valentine, Texas, 

Rarthquake (Appendix IV) 

Deformation possibly related to an earthquake in Texas has been revealed 
from analysis of precise leveling data collected by the National Geodetic 
Survey. Relative elevation changes deduced from first-order leveling surveys 
conducted between 1910 and 1957 show that at least two areas in western 
Texas near the epicenter of the 1931 Valentine earthquake (M = 6.4) subsided 
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relative to their surroundings. Apparent subsidence east of Van Horn, Texas 
is attributed primarily to the combined effects of groundwater withdrawal 
and topography- related survey errors. In contrast, relative subsidence near 
Valentine, Texas, of 11.2 * 1.0 cm extending over a distance of about 20 km 
does not appear to be due to either near-surface effects or leveling errors 
and thus may represent cosesimic deformation of the Valentine earthquake. 

The predominance of subsidence over uplift indicated by the leveling data 
suggests a normal fault mechanism for the earthquake. This is consistent 
with the Cenozoic structure of the region which is characterized by Basin 
and Range morphology. If the observed subsidence was caused by the 1931 
earthquake, the leveling data suggest that the epicenter for this event 
lies considerably closer to Valentine than originally reported. This new 
location is supported by a recent relocation for the Valentine earthquake, 
and by local intensity reports. 

Ill* Detection of Vertical Fault Displacements by Precise Leveling in 

Western Kentucky (Appendix V) 

Relative vertical displacements of bench marks in extreme western 
Kentucky have been determined by comparison of successive leveling surveys 
in 1947 and 1968. The resulting pattern of apparent surface deformation 
shows a steep offset which can be closely modeled by a normal fault buried 
in an elastic half-space. The offset is located near the northern boundary 
of the Mississippi Embayment and the New Madrid seismic zone, an area where 
faults have previously been inferred on the basis of both geological and 
geophysical evidence. If the apparent movement is due to slip along a 
fault, several lines of evidence (regional structure, earthquake data, 
and lineations) suggest that the postulated fault trends north -northeast. 


Thirteen earthquakes were recorded in this area between the times of leveling; 
focal mechanisms exist for three of these. The nearest of these three focal 
mechanisms to the leveling offset implies normal faulting. The magnitude of 
the earthquake, however, appears to be too small to account for the amount of 
slip required by the fault model. Thus the apparent deformation may have 
accmulated with several undetected small earthquakes, or gradually as aseismic 
creep. 

IV. An Empirical Basis for Defi n ing Releveling Anomalies; Implications for 

Crustal Deformation in Southern California 

Standard error estimates for precise leveling are based on loop mis- 
closures and on the repeatability (usually two determinations) of height 
differences between consecutive benchmarks. Standard error estimates as such 
contain little information on systematic errors. Primarily as a result of 
crustal movement studies (e.g. Savage and Church, 1974; Brown, 1978; Citron 
and Brown, 1979; Chi et al., 1980; Jackson et al., 1980; Strange, 1981; 
Reilinger and Brown, 1981) it is now known that systematic leveling errors 
can surpass standard error estimates over distances of only a few kilometers. 
IVhile systematic errors can be recognized in some instances (e.g. Brown, 1978; 
Chi et al., 1980; Ni et al., 1981; Reilinger and Brown, 1981), in many cases 
erroneous measurements cannot be unequivocally identified (e.g. Citron and 
Brown, 1979; Mark et al., 1981). 

Because of existing uncertainties as to the cause or magnitude of sys- 
tematic leveling errors we have begun developing an empirical basis for 
identifying anomalous releveling signals. This project utilizes the stan- 
dardized U.S. releveling data base (all movement profiles plotted on common 
scales) developed at Cornell (Fig. 1). Each profile of apparent elevation 
change (Fig. 2a shows an example from our data base of a standardized ele- 
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vation change profile) is approximated by a series of linear tilt segments which 
are determined by the particular character of the profile (e.g.. Fig. 2b). 

To minimize the influence of unstable benchmarks, each segment must be defined 
by at least S benchmarks and must extend at least 10 km along the route. The 
relative movement between the end points of each segment is then catalogued 
along with the information listed in Table 1 (e.g. location, segment length, 
topography, benchmark spacing, orientation, time interval, etc.). This 
geodynamics data base will allow quantitative determination of the dependence 
of apparent movements on any combination of the parameters listed in Table 1. 
Perhaps most importantly, it will be possible to determine what constitutes 
an artomalous signal both in terms of the magnitude and spatial extent of 
observed movements. 

There are a variety of problems directly relevant to earthquake prediction 
which can be addressed with this new data base. For example, the dependence of 
tilt on topography can be evaulated as a function of tectonic setting (i.e. are 
apparent movements which correlate with terrain confined to tectonically active 
areas such as Southern California or do such correlations exist in presumably 
stable regions?). Addressing this question should help determine to what ex- 
tent correlations with terrain, which clearly occur in Southern California, 
result from real movements or systematic errors. Similarly we will investigate 
the differences and similarities (tilt magnitude and wavelength) between 
apparent movements in contrasting tectonic settings. Such an investigation 
may reveal characteristic "signatures" which are unique to seismically 
active regions and which may prove useful for seismic hazard evaulation. 

Although this project is still in the development stage, we have con- 
ucted a preliminary study incorporating all U.S. profiles of apparent ele- 
vation change in areas of significant relief (slope of ^ 5 m/kra) . Each 
profile was approximated by a series of linear tilt segments. Histograms 
showing the number of segments in each tilt range (.2u rad. range) were con- 
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structed for leveling within the Southern California uplift and for leveling 
throughout the rest of the U.S. (Fig. 3a, b). The following preliminary 
observations can be made from these histograms: 1) in general, tilts within 
the region of the Southern California uplift do not differ substantially 
in magnitude from tilts observed in other areas of significmit relief, 

2) tilts of less than about 2y rad. are difficult to distinguish (at least 
on the basis of tilt magnitude) from "background noise": 3) many of the 
larger tilt anomalies (>2m rad.) have been convincingly associated with tec- 
tonic deformation (e.g., Alaska postseismic deformation, 1952, Kem County 
earthquake deformation, Rio Grande rift, Hebgen Lake, etc.; i.e., this 
approach to characterizing the data base has effectively separated tectonic 
deformation from the "Background noise"); 4) at least some of the Southern 
California releveling measurements fall well above the background level 
and thus may represent tectonic deformation. In regard to this last point, 
it is interesting to note that a number of the releveling measurements de- 
fining preseismic effects of the 1971 San Fernando earthquake (S.F.P.S. on 
Fig. 3a) are characterized by rather anomalous tilts. While some of these 
movements are due to groundwater variations (Rei linger, 1980) others may 
be due to precursory phenomena. Because of their obvious importance to 
earthquake prediction in Southern California, further investigation o^ 
these particular observations is receiving our highe.st priority. 

V. New Releveling Measurements in the Eastern U.S. 

The National Geodetic Survey (NGS) is currently involved in a massive 
releveling program to update the North American vertical datum (Whalen, 1980). 
These new measurements are being incorporated in our crustal movement data 
base as they become available. Figure 4 shows those leveling routes along 
the east coast of the U.S. for which new crustal movement information has 
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been obtained. These recent observations are being used to further investigate 
possible deformation suggested by our earlier studies, Including; uplift of 
the Appalachian Mountains (Citron and Brown ^ 1979); eastward tilting of the 
Atlantic Coastal Plain (Brown and Oliver, 1976); and apparent movements asso- 
ciated with a number of structural features along the east coast (Brown, 1978). 

Perhaps the most interesting result obtained thus far concerns possible 
contemporary deformation in the vicinity of Cape Hatteras, North Carolina. 

Brown (1978) produced and Interpreted profiles of apparent elevation change 
from leveling and tide gauge measurements along the east coast from Calais, 
Maine to Key West, Florida. This study revealed large discrepancies between 
the two techniques which are now believed to be due to systematic errors in 
leveling (Reilinger and Brown, 1981). On the other hand, individual movement 
features along this coastal line correlate with geologic structures, and may 
reflect tectonic deformation (Brown, 1978). Steep tilting is indicated 
near a fault bounding the Connecticut graben, subsidence is noted at the 
Chesapeake embayment. Sharp subsidence near Cape Canaveral and relative up- 
lift near Cape Fear are in regions where geomorphic and seismic evidence 
suggest recent movements. The steepest regional tilting indicated by 
the coastal profile occurs near Cape Hatteras, N.C. Figure 5 shows apparent 
elevation changes along the Cape Hatteras segment of the coastal route 
incorporating the 1979 releveling measurements. As indicated by Figure 5, 
the new survey shows continued tilting in the same sense and at approximately 
the same rate as the earlier measurements. The progressive deformation 
indicated by multiple surveys, the large magnitude of the observed tilting 
and the lack of significant terrain along the route (ruling out elevation 
dependent errors) is strong evidence that these apparent movements represent 
crustal deformation and not measurement errors. The long wavelength of the 


10 


observed tilting argues against near surface effects and for a deep seated tec- 
tonic mechanisn. Our analysis of these new observations, which is currently 
in progress, involves testing the leveling iaeasurcmcmts for internal consis- 
tency and incorporating other geophysical and geological information on neo- 
tectonic activity (e.g. tide gauge measurements, deformation of coastal 
terraces, recent seismic activity, etc.)* Should the measurements stand up to 
this analysis, they will constitute the most persuasive evidence for neo- 
tectonic deformation in the east yet uncovered by leveling. As the east 
coast of the U.S. has experienced major historic earthquakes (e.g. Charleston, 
1886), the possibility that anomalous deformation in the Cape Hatteras area 
reflects a potential seismic hazard will be carefully examined. This example 
illustrates the importance of the new leveling measurements being collected 
by the NGS for evalurt^iq the current state of tectonic activity in the U.S. 
Incorporating these measurements into the existing data base will represent 
a significant part of our future effort. 
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generally the last survey along route) 

rclevcling year - time of releveling (note: if a route has been rcleved more 

than once, each relevcling defines a different segment). 
Southern California code - this code identifies profile segments located within the 

area of the Southern California uplift. 
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Figure 2b: Example of tilt segments used to approximate elevation 

change profiles. 
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(see Figure 4 for location) . 
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NBOTECTONIC DEFORMATION. NEAR-SURFACE 
MOVEMENTS AND SYSTEMATIC ERRORS 
IN U.S. RELBVELING MEASUREMENTS: 
IMPLICATIONS FOR EARTHQUAKE PREDICTION 

by 
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Abstract . Analyses of U.S. releveling 
measurements indicate that derivative crustal 
movement estimates may reflect tectonic 
deformation, near-surface movements, and/or 
systematic errors. Discriminating the 
contributions of these factors is especially 
crucial for unambiguous geodetic detection of 
possible precursory seismic deformations. VDiile 
reliable leveling measurements of co-seismic and 
post-seismic movements are well documented for 
some of the larger (M > 6) dip-slip earthquakes, 
leveling evidence for pre-seismic motion is 
generally sparse and often ambiguous. Subtle 
earthquake-related motions may be masked by both 
aseismic movements and systematic errors. Deep 
magma injection and surficial groundwater 
withdrawal are two mechanisms which are shown to 
cause surface movements which, under some 
circumstances, could be misidentif led as seismic- 
related. Of more concern, perhaps, are 
systematic measurement errors. Topography 
dependent errors are an exceptionally troublesome 
type, perhaps affecting as much as 20% of U.S. 
leveling. However, other varieties of systematic 
error also contribute to the uncertainty. 
Discrepancies between leveling and tide gauge 
data and within nets of leveling alone suggest 
large, long baseline accumulations of error. In 
many cases, aseismic and erroneous contributions 
can not be unequivocally determined _ex pos t 
facto . However, a comprehensive examination of 
the NGS crustal movement data base, representing 
a large sampling of the entire U.S. Level Net, 
provides perspective and criteria needed to begin 
to recognize movement directly related to 
earthquake activity. 

Perhaps the most extensive set of measurements 
relevant to earthquake activity exists in 
southern California, where much attention has 
recently been focused. Reevaluation of some of 
these leveling observations indicates that while 
some appear to reflect tectonic deformation, 


others are suspect because of indications of 
systematic errors and near-surface, non-eectonic 
movements. Specifically, possible preselsmic 
movements reported for the 1971 San Fernando 
earthquake in the vicinity of the earthquake 
fault as well as approximately 30 km northwest of 
the epicenter may bo due to systematic errors. 
Movements near the San Gabriel fault, initially 
ascribed to the Palmdale Bulge and more recently 
to preselsmic effects of the 1971 San Fernando 
earthquake apparently reflect near-surface 
sediment compaction due to water table 
fluctuations. Similarly, there Is strong 
evidence of contamination by rod calibration 
errors la some of the releveling observations 
used to define the southern portion of the 
Palmdale Bulge (Llano to Azusa, California). The 
reality of the Palmdale Bulge itself must be 
questioned in view of this reevaluation. In 
contrast, possible tilting southwest of Palmdale 
between 1961 and 1964 is not easily related to 
systematic errors or near-surface movements and 
thus may represent tectonic deformation. Whether 
this tilt anomaly was due to preselsmic effects 
of the San Fernando earthquake or a mechanically 
separate tectonic event Is presently unknown. 


Introduction 


Vertical movements of the earth's crust are 
commonly expected to accompany the various phases 
of strain buildup and release associated with 
major earthquakes. Observations of vertical co- 
selsmic and post-seismic movements using precise 
lovol Ing nri* wi-l l-<locnmc»ntoil. Ilowovor, repr)rtH 
ot prusulsiulc muvumunt in Che U.S. are rare uiiU , 
as will be argued, questionable. Recognition of 
true pre-scismlc motJ is complicated by 
systematic leveling errors, near-surface non- 
tectonlc processes (e.g., fluid withdrawal), the 
general lack of sufficiently redundant and 
extensive surveys, and the fact that significant 
changes in elevation have been identified which 
are unrelated to earthquakes. Such "noise" could 
easily hide a pre-seismlc signal. Considerable 
uncertainiv exists as to the extent and magnitude 
of these obscuring "movements". Direct 
determination of their effects is often extremely 
difficult. However, some perspective on these 
problems can be obtained froi? empirical analyses 
of existing National Ceodetlc Survey (NGS) 
relevellng (Figure 1). Such an analysis forms 
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Che basis of chLs reporc. 

This paper reviews some factors Chat muse be 
considered when atcempuing Co extract tectonic 
j.:>formation, especially chose relevant to 
earthquake predicCion, from historic releveling 
observations. ivvidence on the extent and nature 
of systematic, errors, non-Cectonic movements, and 
tectonic deformation (both earthquake related 
deformation an<1 tectonic movements unassociaCed 
with earthquakes) .'Jvom U»S. releveling 
measurements is presented > Specific criteria to 
help recogni/.c suspect movements are developed 
and illustrated by application Co a reevaluation 
of '•ertain southern California leveling results 
of parciculnr interest in earthquake prediction* 


Systematic Errors in Leveling 


Ac Che root of much of Che current debate 
regarding level ing-derived estimates of ustal 
motion is the prevailing i ncertalnCy as o Che 
role of systematic measurement errors. In 
particular, systematic errors which accumulate 
with relief have become a central Issue In 
crustal movement research . While errors of this 
type have been known to geodesists for some time 
(e.g., Bomford, 1971), their Influence has been 
considered too small to be of concern in most 
geodetic applications. However, new field 
experiments carried out by the NGS (Whalen, 1980) 
and empirical analyses (Brewn et al., 1980) 
confirm earlier suspicions (e.g.. Savage and 
Church, 1974; Brown and Oliver, 1976; Citron and 
Brown, 1979; Jackson et al., 1980; Chi et al., 
1980) that topography-induced systematic errors 
are larger and more common than heretofore 
established and consequently that such errors can 
be and probably have been misinterpreted as 
tectonic motions of the crust. 

Topography-correlated errors can arise from 
Improperly calibrated leveling rods and from 
unequal atmospheric refraction of the foresight 
and backsight readings* The effects of these two 
sources of error should differ in a number of 
respects and thus In principle may be 
distinguished. For example, ficticious movements 
resulting from rod calibration errors should 
correlate rather closely with detailed 
topography, and change magnitude only where rod 
pairs are changed. In contrast, atmospheric 
refraction will be Independent of the rods used 
in the survey. Furthermore, refraction errors 
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cnn bu cxpucCinl tu accumulate in a more complex 
manner because refraction depends on a variety of 
parameters which may vary significantly during 
the course of a given survey (e.g., near surface 
temperature gradients, individual sight lengths, 
wind, etc.)* Because of procedural changes (a 
tendency towards shorter sight lengths In newer 
surveys), atmospheric refraction should more 
often than not result In ficticious moveuents 
which show a positive correlation with topography 
(l.e., high areas will appear to be rising) while 
errors due to rod mlscallbratlon should have no 
preference for positive or negative correlations. 
In practice, these distinctions are not always 
easy to draw. However, a preliminary survey of 
NGS relevel Ing estimates of elevation change 
vhlch correlate with topography Indicates that 
about 7SZ display positive correlations. VAille 
rod calibration errors appear to effect some 
leveling observations (e.g. Jackson et al., 
1980; this paper) this preliminary result 
suggests that refraction may be the more 
pervasive source of elevation correlated error. 

The expected magnitude for refraction error is 
a point of considerable uncertainty. According 
to one approximation (l.e., Kukkamakl , 1938), 
this error is proportional to the height 
difference between benchmarks, the temperature 
difference between 0.5 m and 2.5 m above the 
ground, and the square of the sight length used 
In the observation. Figure 2 shows 
representative values for refraction error using 
constants given by Uoldahl (1980). For 
reasonable temperature differences (1-2 degrees 
C: Whalen, I960) and sight lengths (25-75 m) , 
errors as large as 30-40 mm or more can easily 
accumulate over height differences of 100 ro 
(300-400 ppm). Since refraction error will 
usually have the same sign. Its effect should 
tend to cancel when differencing surveys to 
compute moveiiieiii; . This rationale has often be»*n 
used as an urgument for Ignoring the effect. 
However, If surveys are conducted using different 
sight lengths (Strange, 1981) and/or under 
different mic rometeorologlcal conditions, the 
refraction effect will result In what appear to 
be movements that roughly correlate with relief. 

Examples of apparent movement correlating with 
elevation are numerous (Brown et al., 1980). 
Approximately 20Z of U.S. relevellng observations 
show visual correlations between apparent 
movement and topography. The magnitude of the 
effect often reaches 30-40 mm per 100 m change In 
elevation. Although, correlation with topography 
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oloitc 1 h InMiifC Ic Lent to warrant rejection of a 
tectonic interpretation (e«g«, see Reilinger et 
al., 1977), it is clearly grounds for suspicion 
particularly when correlations persist at short 
wavelengths (e.g.. Figures 18 and 19) or when 
multiple surveys give Inconsistent results. For 
example, Figure 3a shows apparent vertical 
movements and terrain along the route from 
Colorado Springs to Leadville, Colorado based on 
surveys conducted in 1925, 1953, and 1954. The 
reversal of the 1953-1925 apparent tilt during 
the period 1954-1953 strongly suggests systematic 
error, and the correlation with terrain (Figure 
3b) points toward an elevation correlated error. 
Since elevation-dependent errors may contaminate 
a significant portion of the NCS releveling data 
base, the possibility of such errors must always 
be considered prior to invoking tectonic 
explanations for apparent movements which 
correlate with relief. 

Topography dependent error is not the only type 
of systematic error affecting U.S. leveling 
measurements. Table 1 lists a number of areas 
where comparison of repeated leveling 
measurements show large systematic discrepancies 
which may be due to errors in the observations. 
These examples occur in areas of generally 
subdued relief, thus ruling out elevation- 
correlated errors. If leveling errors are solely 
responsible for these discrepancies, whatever 
their cause, they range in magnitude from .3 
mm/km to 1 mm/km and remain systematic (i.e., 
accumulate monotonically) for distances ranging 
from about 100 km to over 600 km. 

Figure 4 shows three different estimates of 
elevation change (assuming constant rates of 
movement over the time period between levelings) 
along the east coast of the U.S. from Maine to 
Florida: 1) from unadjusted leveling 
measurements; 2) from tide gauge records 
(squares); and 3) from the same leveling 
observations adjusted with standard least squares 
procedures for consistency with other repeated 
leveling lines which form circuits extending 
inland from tliu coast (tide gauge data were not 
used in the adjustment, Jurkowski et al., 1979). 
The leveling measurements span an approximately 
30 yr. time interval. The fact that the relative 
movement between Maine and Florida is 
substantially reduced through the adjustment and 
the fact that the adjusted leveling profile is 
more consistent with tide gauge data (although 
serious discrepancies still remain) indicate that 
the regional north-south tilt results from 
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systcmaClc crrurti in Che leveling observations. 
The error rumains more or less sysCcmaCic over 
discances of lOOO's of kilometers, and on some 
sections (e.g., 1300-2600 km) reaches .5 mro/km. 

Comparison of leveling and tide gauge estimates 
of crustal movement along Che west coast of the 
U.S. between Astoria, Oregon, and Crescent City, 
California show similar discrepancies (Brown et 
al., 1980). Unlike the east coast profile, 
apparent crustal movements along Che west coast 
were derived from only two surveys and were thus 
not subject to possible temporal bias due Co 
striwging together segments covering different 
time intervals. For Che west coast profile, the 
north-south error reaches .3 mm/km and remains 
systematic over a distance of 380 km. 

The cause of the apparent errors in these 
coastal surveys is presently unknown. The 
predominantly north-south orientation of the 
coastal profiles may suggest unequal lighting or 
ocher factors which are believed Co 
preferentially accumulate on north-south lines 
(Bomford, 1971). However, the substantial 
reduction of Che apparent tilt indicated by the 
east coast profile when adjusted with inland data 
suggests that the error may be related Co the 
proximity of Che leveling route Co the coast 
(i.c., the error did not effect, nr had less of 
an effucC, on pr4>filcs further inland). 

Suspect movements are not restricted to coastal 
ar.d north-south profiles. For example, consider 
Che large apparent tilt across Che U.S. 
midcontinenC identified by Brown and Oliver 
(1976) from relcveling between Davis Junction, 
Illinois and Willard, Ohio (Figure 5). This tilt 
is perhaps the largest apparent movement defined 
by leveling in the eastern U.S. The tilt anomaly 
shows no clear relationship Co geologic structure 
and is inconsistent with movements inferred from 
comparisons of water level gauges In the great 
lakes (Brown and Oliver, 1976). Figure 5 shows 
Che results of a loop closure analysis (see Chi 
et al., 1980 for discussion of method) for 
circuits including the Davis Junction Co Willard 
route. The fact that misclosures are 
considerably larger when the circuits are closed 
with Che I9f)/-I909 surveys between Dnvls Junction 
and Willard than with the 1930-1947 surveys, even 
Chough the remainder of the loop was surveyed 
more closely in time to the 1967-1969 Interval, 
suggests Chat the large apparent tilt of the 
interior plains is due Co systematic error and 
not real ground motion. Such an error would have 
CO resell I mm/kiQ and remain systematic for 
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discanccH of wull over 500 km. 

la both chu coastal and interior examples cited 
above, the discrepancies are characterized by 
consistent accumulation over large distances. 
The effects of these long baseline discrepancies 
are often reduced below the level of concern for 
most geodetic applications by network 
adjustments. However, when examined in the 
context of evaluating crustal movements such 
discrepancie.' become quite significant. Although 
the net apr irent movement over a profile can be 
large, the tilt rates resulting from these errors 
are rather low, especially when compared with 
those exhibited by unequivocal examples of real 
movement. Tilt may therefore be the more 
diagnostic parameter in evaluating reliability of 
crustal movement estimates. 


Near-Surface Movements 


In addition to systematic errors, releveling 
measurements are influenced by near-surface 
movements which can mask, or be mistaken for 
deep-seated tectonic motion. Near-surface 
effects include: a) Benchmark instability (frost 
heave, soil moisture and temperature changes, 
human disturbance), b) Surface failure (land 
slides, mine and cavern collapse), c) Loading 
(reservoir Impoundincnt , building settlement), and 
d) Fluid withdrawal (water, oil, gas). Benchmark 
instability and surface failure are often easily 
identified or are of such local extent that they 
are not a serious problem for regional tectonic 
studies. Such effects can, however, complicate 
local investigations - for example, of movements 
near earthquake faults. Ncar-surface soil or 
sediment compaction due to earthquake ground- 
shaking may be responsible for the predominance 
of subsidence over uplift near many earthquake 
faults (Savage and Hastie, 1966). Subsidence due 
to surface loading and fluid wlthd* al is, in 
general, easily related to human activity. In 
fact, leveling has proven quite effective at 
monitoring such movements, with important 
engineering applications (e.g., Poland and Davis, 
1969). However, near-surface movements, and in 
particular muvumunts due to variations of water 
levels in aquifers, appear to be more widespread 
than previously reported. In addition, such 
effects can be subtle and subsequently 
misidentif ied as tectonic deformation. 

Figure 6 shews NGS releveling profiles in the 
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U.S. whlcVt Imllcnte subsidence relative to 
surrounding areas and which overlie aquifer 
systems which have experienced variations in 
water levels due either to pumping or natural 
causes. These movements may therefore represent 
sediment compaction associated with these water 
level variations, and are consequently suspect as 
indicators of tectonic motion. It is Interesting 
to note from Figure 6 that these apparently near- 
surface movements are quite common In southern 
California, affecting much of the area peripheral 
to the Palmdale Bulge. Arguments will be 
presented later suggesting that In at least one 
case such groundwater subsidence In southern 
California has been mlsldentlf led as tectonic 
uplift of the adjacent areas* 

A particular example, not previously reported, 
which Illustrates criteria which can be used to 
recognize near-surface sediment compaction Is the 
relative subsidence of the Los Angeles basin. 
Figure 7 gives a map of the L.A. Basin showing 
contours of basement depth and the location of a 
leveling route traversing the basin. Also shown 
are elevation changes along the leveling route 
and the history of water level decline measured 
In an observation well near the center of the 
basin. Subsidence near the center of the basin 
reaches IS cm relative to the periphery and 
extends over i distance of 40 km. The observed 
subsidence correlates spatially with aquifer 
geometry and temporally with the history of water 
level decline. In addition, Che magnitude of Che 
effect (i.e., Che ratio of subsidence to water 
level decline) Is comparable to observations In 
other areas (Poland and Davis, 1969). Had the 
relationship between aquifer geometry and 
subsidence not been noticed. It Is possible Chat 
these measurements could have been misinterpreted 
as tectonic motion. 


Vertical Movements and Earthquakes 


In spite of the substantial difficulties 
associated with relevellng estimates of crustal 
movement, some of which have been described In 
previous sections, the capability of Che leveling 
technique for monitoring tectonic earth movements 
Is well established. In a number of cases, 
relatively subtle earth movements (I.e., tilts 
few X 10”O rad and tilt rates few x ICT® rad/yr) 
have been identified. In this section we briefly 
review relevellng evidence for earthquake related 
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deformation In the U.S. and use specific examples 
to Illustrate some of the criteria employed to 
Identify real tectonic movements. 

The best examples of tectonic deformations 
measured by leveling In the U.S. are those In 
the vicinity uf major earthquakes. Figure 1 and 
Table 2 review those U.S. earthquakes for which 
vertical movements have been reported. All of 
these earthquakes are associated vrlth faults that 
have a significant component of dlp-sllp movement 
(with the possible exception of the 1940 Imperial 
Valley earthquake). Up to the present, there Is 
no clear evidence from U.S. rcleveling 
measurements for permanent vertical deformation 
associated with purely strike-slip faulting 
although present observations are not sufficient 
to rule such out . 

The most obvious vertical movements are those 
accompanying the earthquake (coselsmlc). 
Coselsmlc deformation has been well-documented 
for several of the larger (M >6) normal and 
thrust earthquakes which have occurred In areas 
of preexisting geodetic control (Table 2). 
Observed movements range In magnitude from a few 
cm to a few m depending on the size of the 
earthquake and the proximity of the leveling 
measurements to the eplcentral area. In general, 
coselsmlc movements are well explained by elastic 
dislocation theory (Savage and Hastle, 1966) 
although complications can arise from such 
factors as near-surface soil or sediment 
compaction due to ground shaking. 

Post-selsmlc vertical movements have also been 
observed by rulcvellng for some of the larger 
dlp-sllp earthquakes (see Table 2). These 
movements are usually smaller than associated 
coselsmlc movements; however like coselsmlc 
movements they can often be Identified by their 
close spatial and temporal association with 
earthquakes, and In some cases surface faulting. 
Where sufficient observations exist, post-selsmlc 
deformation rates appear to decrease 
exponentially from the time of the earthquake. 
For example, movements near Anchorage following 
the 1964 Alaska earthquake are shown In Figure 8 
(Brown et al., 1977). The Alaska earthquake, one 
of the largest events ever recorded, occurred 
where the oceanic Pacific plate is being thrust 
under the continental North American plate at a 
rate of over 5 cm/yr (Plafker, 1972). Savage and 
Hastie (1966) and Hastle and Savage (1970), using 
a dislocation model of thrust faulting, showed 
that the coselsmlc displacements were consistent 
with low-augli- thrusting. Fust-seismic movements 
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near Anchorage (Figure 8) amounted to as much as 
0.55 m of land uplift at an exponentially 
decreasing rate during the decade following the 
earthquake. Additional evidence for deformation 
following the Alaska earthquake was reported by 
Prescott and Lisowski (1977) from analysis of 
detailed leveling arrays on Middleton Island in 
the Gulf of Alaska. Tilts associated with the 
Alaska post-seismic movements were v^n the order 
of 10“^ to 10~” rad. There is still considerable 
debate as to the mechanism responsible for post- 
seismic movements t but at least some of the 
observations appear consistent with after-slip on 
the fault, or an extension of the fault that 
ruptured during the earthquake, although other 
explanations have been proposed (e.g., Nur and 
Mavko, 1974; Wahr and Wyss, 1980). 

While co-scismic and post-seismic movements are 
well established for at least some earthquakes, 
clear evidence for preseismic deformation from 
U.S. relevcllng measurements is quite rare. 
This may be due to a lack of appropriate 
measurements as opposed to the absence of such 
movements since it is unusual to have multiple 
levelings of sufficient proximity prior to an 
earthquake. Precursory vertical movements have 
been suggesteu from leveling measurements for 
only three U.S. earthquakes: the 1959 magnitude 
7.1 Hebgen Lake, Montana, the 1971 magnitude 6.4 
San Fernando, California and the 1973 laagnitude 
6.0 Point Mugn, California earthquakes (see Table 
2). The evidence for preseismic movement near 
Point Mugu is marginal, both because the 
proposed movements are barely significant 
relative to random error estimates and because 
the area was subject to surficlal subsidence due 
to groundwater withdrawal during the period of 
Interest (Cnitlli* et al., 1977). The 1971 San 
Fernando earthquake is exceptional in that 
significant releveling was available for the 
eplcentral area prior to the earthquake. These 
observations were analyzed after the earthquake 
and were Interpreted to indicate precursory 
movements (Castle et al., 1974). However, 
reevaluation of the relevant leveling 
observations, described in a later section of 
this paper, cast some doubts on the reliability 
of these measurements and hence on their tectonic 
significance. Reillnger et al. (1977) found 
evidence for possible precursory uplift 
throughout a broad region surrounding the area of 
major co-seismic movement of the 1959 Hebgen Lake 
earthquake which apparently accumulated at a rate 
of 3-5 mm/yr (Figure 9). The zone of uplift is 
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defined by five Independent elevation change 
profiles derived from 12 independent surveys. 
Although three of the five movement profiles show 
positive correlation with topography (i.e. high 
areas going up), one shows a negative correlation 
and one shows no correlation; yet all indicate a 
consistent sense of movement. This consistency 
argues strongly against elevation correlated 
errors as the cause of the observed uplift. The 
doming stands out distinctly In relation to 
movements In surrounding areas, and shows a close 
spatial correlation with the zone of major co- 
selsmlc deformation and aftershock activity for 
the 1959 earthquake (Figure 9). In addition, the 
geodetically measured deformation is consistent 
in sign with Cenozolc deformation deduced from 
geologic structure (Relllnger et al., 1977). 
Tilts associated with this uplift range from 3 - 
7 X 10~^ rad with associated tilt rates between 1 
-3 X 10"^ rad/yr. Although Relllnger et al. 
(1977) suggest that doming began prior to the 
earthquake, because of the limited number of pre- 
earthquake leveling measurements, it is 
impossible to prove that the activity was 
precursory (i.e., doming may have accompanied, 
and/or immediately followed the earthquake). In 
sum, therefore, leveling evidence for vertical 
movements preceding any U.S. earthquake Is weak 
in both quantity and quality. 


Other Tectonic Deformation 


Recognizing real tectonic deformation from 
relevellng, although necessary, is not sufficient 
grounds to Infer that they are directly relevant 
to the cartlu|uake prediction problem. Earthquake 
related movements must be separated from 
movements due to other deep seated processes, 
such as Inelastic deformation (fault creep, 
folding; e.g. Thatcher, 1931), isostatic 
adjustments and magmatic activity. These 
mechanisms are believed, on the basis of 
observational evidence, to result in contemporary 
vertical movements which are sufficiently rapid 
to be detected by releveling measurements. 

Movements due to subsurface magmatic activity, 
for example, are not restricted to volcanically 
active regions (e.g., Hawaii, Iceland, Japan), 
having been reported in Yellowstone National Park 
(Relllnger et al., 1977; Pelton and Smith, 1979), 
and the Rio Grande rift (Relllnger et al., 1980) 
as well. Crustal uplift in the Central Rio 
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Grande rife ilLustracea cecConlc deformation 
which appearH to be unrelated to major earthquake 
activity* Thu existence of an active magma body 
beneath the central Rio Grande rift was inferred 
primarily on the basis of geophysical, and some 
geological information (Sanford et al*, 1977). 
The magma body is believed to consist of a thin 
sill at a depth of about 20 km (Figure 10). 
Elevation change profiles along the routes shown 
in Fig<;re 10 are given in Figure 11. All three 
profiles indicate uplift of the area overlying 
the magma body. The observed uplift is believed 
to be due to tectonic deformation and not 
measurement errors or near-surface movements 
because: 1) uplift is defined by three 

independent elevation change profiles; 2) while 
the two east-west profiles show a rough negative 
correlation with topography near the area of 
uplift, the north-south profile shows no 
correlation, thus ruling out elevation-dependent 
errors as thu primary cause of the observed 
movements; T) the Helen to Amarillo profile 

demonstrates that the uplift of the rift is 
anomalous relative to points to the east; 4) 
georaorphic evidence for post-Pliocene deformation 
(Bachman and Mehnart, 1978) is consistent in sign 
with the geodetic observations; 5) anomalous 
uplift occurs directly above the magma body; and 
6) modeling studies indicate that uplift could 
result from activity within the magma body. If 
uplift is accumulating more or less continuously 
as suggested (Reilinger et al., 1980), it is 
characterized by an average rate of 4 mm/yr^ 
with corresponding tilt rates of 5 to 10 x 10~ 
rad/yr. In spite of the rather compelling 
independent evidence for an active magma body 
beneath the area of uplift, it has been suggested 
that these movements may be due to an impending 
earthquake (Koseluk and Bischke, 1981). 
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Reliability Criteria 


The selected cases described above demonstrate 
both the utility and limitations of geodetic 
leveling to detect tilts of a few x IGT” rad and 
tilt rates of a few x 10~^ rad/yr. Thus while 
non-tectonic influences (e.g« systematic error) 
can obscure real earth movement, the technique 
has clearly proven effective at monitoring 
relatively subtle tectonic deformation* It Is 
essential, however, that individual releveling 
observations be exami''.ed in detail for possible 
contamination by systematic errors and near 
surface movements prior to invoking tectonic 
explanations. Particularly effective 
quantitative techniques include comparison of 
forward and reverse levelings (e.g.. Savage and 
Church, 1974) and loop closure analysis (e.g., 
Chi et al., 1980). In addition, the following 
qualitative criteria, some of which were 
illustrated by the previous examples, have proven 
useful for evaluating the reliability of 
particular data sets (Brown et al*, 1980): 1) 
magnitude of apparent movements relative to 
possible errors (since many errors remain poorly 
understood this is equivalent to determining 
whether the movements in question stand out in 
relation to "background noise"); 2) consistent 
temporal behavior when multiple levelings are 
available (e.g., Alaska); 3) relations with 
independent geophysical or geologic estimates of 
recent movement (e.g., tide guage, lake levels, 
tilt meters, horizontal movements, geomorphic 
evidence, etc.); 4) consistent movements when 
multiple leveling lines cross a given feature 
(e.g., Hebgen Lake, Rio Grande rift); 5) 
correlation with geologic structure and tectonic 
activity (e.g., Hebgen Lake ); 6) lack of 
correlation with topography ruling out possible 
elevation correlated errors; 7) lack of 
relationship to possible near-surface processes 
(e.g., fluid withdrawal, reservoir impoundment, 
etc.); 8) lack of relationship between apparent 
movements ami procedural changes (changes In 
sight lengths, rod or instrument changes); and 9) 
consistency of inferred mechanism with tectonic 
setting (e.g., Alaska). 

Although no single criteria is sufficient to 
warrant acceptance or rejection of a set of 
observations, when multiple criteria give 
consistect results they provide strong evidence 
for the reliability of the measurements in 
question* 
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A Case Study: Southern California Relevellng 
Meaeuremcnte 


Much attention has recently been focused on 
leveling in southern California, where there is 
both considerable concern about future 
earthquakes and an abundance of leveling 
observations. Using the above reliability 
criteria, developed through analysis of the much 
broader data base of U.S. relevellng, we have 
reevaluated some of the observations used to 
deduce presclsmlc movements for the 1971 San 
Fernando earthquake as well as the Palmdale 
Bulge. Our reevaluation, representing a 
different perspective, suggests that many of the 
southern California measurements are 
significantly affected by both topography- 
dependent errors and near-surface movements. On 
the other hand, at least some of the observations 
may still reflect tectonic deformation. Thus, 
the configuration of the Palmdale Bulge will, at 
the very least, require revision in light of 
improved understanding of those factors which can 
influence relevellng measurements. 

In our analysis of southern California 
relevellng obuerv.itlons, data have been displayed 
in terms of relative movements for sequential 
time intervals along the pertinent segments of 
the leveling routes. This contrasts with 
previous attempts to tie the observations to a 
tide gauge In order to relate movements to sea 
level. Analyzing relative movements minimizes 
the effects of systematic errors, which can 
accumulate to rather substantial amounts over the 
100-200 km distance to the tide gauge, and as 
will be demonstrated, greatly simplifies 
interpretation of the observations. 

Figure 12 shows those leveling routes in 
Southern California for which crustal movement 
information has been investigated for this study. 
This information was used to define both the 
configuration of the Palmdale Bulge and 
preseismlc movements for the 1971 San Fernando 
earthquake. Possible preseismlc movements were 
reported by Castle et al. (1974) in the vicinity 
of the eartliquake fault (segment I) and 30 km 
northwest of the epicenter (segment II). Strange 
(1981) subsequently reported evidence of 
preseismlc movements Just north and east of 
Saugus (along segment II and III). These 
observations were previously used to define the 
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Palmdale Bulge (Cascle et al*, 1976; Castle, 
1973). Apparent movements along segment 111 and 
along the Llano to Azusa route (Figure 12f) were 
also instrumental in determining the 
configuration of the Palmdale Bulge (Castle et 
al., 1976; Castle, 1978). Each of these features 
will he discussed separately below. 

The sequence of movements along segment 1 
crossing near the area of surface faulting are 
shown in Figure 13. Coseismic movement 
consisting of relative subsidence south of the 
San Fernando Fault and relative uplift north of 
the fault arc clearly indicated for the 1969-1971 
interval. These movements are roughly consistent 
with elastic rebound accompanying thrust faulting 
(Savage et al., 1975). Possible preseismic 
tilting up to the north is indicated by the 
profiles for the time intervals 1955-1961, 
1961-1964, and 1964-1965. Apparently no tilt 
accumulated along this section betwoen 1965-1969. 
Figure 14a shows relative movements between 
points near the ends of this profile segment 
plotted as a function of time. The temporal 
consistency of these movements is, in itself, 
normally evidence that the measurements reflect 
real movements (criteria 2). However, there are 
two reasons to suspect systematic error, and in 
particular refraction errors, rather than true 
ground motion. 

Examination of the profiles in Figure 13 
indicates that the observed tilting correlates 
with topography (criteria 6). This correlation, 
although suggestive, is not sufficient to confirm 
systematic error because real movements can in 
some cases correlate with relief (e.g., Reilinger 
et al., 1977). However, the sequence of apparent 
tilts between 1955 and 1969 show a systematic 
relationship to the sight lengths used for 
different surveys (Figure 14b); a relationship 
that is consistent with that expected from 
refraction errors (Holdahl, 198C) (criteria 8). 
The 3°C temperature difference that results in a 
good fit to the observations, although somewhat 
higher than dally average temperature differences 
observed in other areas (1 -2°C, Whalen, 1980), 
may oot be unreasonable for the spring and summer 
months in southern California (Holdahl, 1980 
reports that temperature differences "may 
frequently aCLaln values up to 4°C shortly after 
noon during the summer"). Although Strange 
(1981) contends that significant tilting persists 
along this segment of the profile in spite of 
refraction corrections, the effectiveness of 
corrections made without on site temperature 


mean II cements Is not well estnbllshed. In view of 
Che possibility of refrscClon errors* the 
tectonic significance of the sequence of apparent 
tilts shown In Figure 13 remains ambiguous. 

Figure IS shows profiles of relative elevation 
change for sequential time Intervals crossing Che 
two areas of reported preselsmic deformation 
northwest of the epicenter (Segment ll; Figure 
12). During the 1953-1964 Interval the main 
deformation consisted of subsidence In the 
vicinity of Saugus relative to points farther 
north (ruled area of plot). This movement was 
originally attributed to the Palmdale Bulge 
(Castle et al., 1976; Castle* 1978) and more 
recently to preselsmic effects of the San 
Fernando earthquake (Strange* 1981). However* 
analysis of relevellng measurements throughout 
the Saugus Basin Indicates that relative 
subsidence shows a close correlation with the 
geometry of the Saugus aquifer and the history of 
water level decline (Bellinger* 1980). The 
spatial correlation between the zone of 
subsidence and the Saugus aquifer is shown in 
Figure 16. In addition to the spatial and 
temporal correlation* the degree of subsidence of 
individual benchmarks %d.thin the aquifer is 
roughly proportional to the product of aquifer 
thickness and water level decline (Figure 17), 
the expected relationship for near surface 
compaction (e.g. Terzaghl and Peck, 1967). These 
observations strongly suggest that subsidence 
above the Saugus aquifer results from near- 
surface sediment compaction due to fluctuations 
of the water level and not from tectonic 
deformation (criteria 7). This result is 
particularly Important to the current controversy 
surrounding the Palmdale Bulge (e.g. Stein, 1981) 
since* unlike many of the measurements defining 
the Bulge* those in the Saugus area do not 
correlate with relief. 

The other large possible movements shown in 
Figure 15 occurred between 1965 and 1968 and 
between T968 and 1969. These observations were 
the primary evidence used to Infer pre-earthquake 
slip at depth on the earthquake fault (Thatcher, 
1976). The 1965-1968 movements consisted of 
uplift of the north section relative to the south 
by about 6 cm. The 1968-1969 movements were* in 
essence* a reversal of the 1965-1968 movements 
(criteria 2). The important point Is that both 
sets of apparent movements were dependent upon 
the 1968 survey. This Is Illustrated by the 
bottom-most plot In Figure 15* which shows the 
general absence of movement for the 1965-1969 
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interval* Therefore, either we were fortunate 
enough to catch preeoisuic deformation at a time 
of significant deflection (1968), and again when 
the movements had exactly reversed themselves 
(1969), or the 1968 survey was in error. Wijile 
oscillatory movements may have occurred, the 
possiMlity of errors in the 1968 survey is at 
least as likely, particularly in light of the now 
suspect results south of the epicenter, and 
similarly suspect trends identified in leveling 
observations In other parts of the country* 

The possibility that refraction errors 
contaminate leveling measurements south of the 
epicenter naturally raises the question as to 
whether this same effect is responsible for the 
apparent error in the 1968 survey northwest of 
the earthquake* The steep tilts Indicated by the 
1965-1968 and the 1968-1969 movement profiles 
occur where there is a corresponding steep slope 
in topography (25 to 40 km)* However, the 
topographic slope is so steep (>*04 rad) that 
only short sight-lengths could be used, making it 
unlikely that atmospheric refraction coupled to 
sight-lengths was a significant effect* Unusual 
near-surface temperature differences at the time 
of the survey, or other elevation-correlated 
errors, such as miscallbrated leveling rods 
(Jackson et al., 1980) may have affected these 
observations (see also Stein, 1981)* 

Apparent movements along the route from Llano 
to Azusa, California were used to define the 
configuration of the Palmdale Bulge (Fig* 12)* 
Elevation change profiles along this route 
derived from surveys in 1934, 1962, and 1971 are 
shown in Figure 18* The 1962-1971 profile 
indie . ' about 10 cm uplift at Llano relative to 
Azus . lowevr.r, the following observations 
stro ./^> vggett. that this apparent movement is 
due ij systematic error, and in particular rod 
calibration error in the 1962 survey (note that 
these are different data than those reported by 
Jackson et al*, 1980): 1) elevation changes 
along the segment of the 1962-1971 profile where 
most of the relative movement occurs (between A 
and B, Fig. 18) show a detailed correlation with 
topography along the leveling route* The 
relationship between apparent movement and 
terrain, which has a correlation coefficient of 
*98 is Illustrated in Figure 19* Such a 
correlation with detailed topographic features is 
suggestive of systematic errors (criteria 6); 2) 
point A in Figure 18 where the correlation 
between raovemont and terrain ends rather abruptly 
is precisuly the point where the leveling 
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instrument nnd rods were changed In the 1962 
survey (criteria 3); 3) the elevation change 
profile derived from the 1934 and 1962 surveys 
along this route indicates movements in the 
opposite sense and of Che same magnitude as the 
1962-1971 profile. This Is clearly Illustrated 
by the 1934-1971 profile (Fig. 18} which shows no 
significant movement. This reversal in Che sense 
of apparent movement casts serious doubts on the 
1962 survey (criteria 2). These three points 
strongly suggest that apparent movements between 
Llano and previously used to define Che 
Palmdale Bulge, actually result from rod 
calibration errors In the 1962 survey. 

Movements along leveling route III in Figure 12 
were also important for determining the geometry 
of the Palmdale Bulge (Castle et al., 1976; 
Castle, 1978). In contrast to the previous 
examples, these apparent movements do not seem Co 
be due Co either systematic errors or near- 
surface effects and thus may represent tectonic 
deformation. 

Figure 20 shows Che sequence of relative 
movements along Che survey route south of 
Palmdale. The major tilt event occurred between 
1961 and 1964 and amounted Co more than 10 cm of 
relative movement over a distance of 20 km. This 
corresponds to a tilt of 5 X 10~ rad. The 
general absence of movements for the 1955-1961 
interval, the 1964-1965 Interval and the 
1965-1971 interval attests to the reliability of 
all of these surveys (i.e., comparison of the 
1955 or 1961 surveys with any of the later 
surveys will give roughly the same result). This 
implies chat Che 1961 to 1964 tilt event is in 
fact defined by five independent surveys. This 
is illustrated by the bottom most profile in Fig. 
20 wlilcli slgnlficuuC tilting for the period 
1961-1965 which, although reduced in amplitude, 
is quite similar to Chat shown by the 1961-1964 
profile (criteria 2). In addition, this tilt 
anomaly does not show a strong correlation with 
topography (i.e., the direction of tilting does 
not reverse where the topographic slope reverses) 
(criteria 6). Furthermore, the sequence of 
relative movements show no relationship either to 
changes in leveling rods or to changes in sight 
lengths (criteria 8). This evidence suggests 
that Che apparent tilting south of Palmdale 
reflects real crustal movements. Whether the 
tilt anomaly was a precur.'^or Co the San Fernando 
earthquake or represented a mechan.lealxy separate 
event is presently unknown. 




PACK 19 


Discussion and Conclusions 


In spice of laudable progress In developing 
sophlsClcaCed new geodetic methods (e.g., VLBI, 
Laser Ranging, GPS) relevellng measurements 
continue to be the most accurate (over 
appropriate distances) and widespread source of 
Information on contemporary vertical movements of 
the continental crust* As such they constitute 
an Important input to the earthquake prediction 
problem* Previous Investigations, a few of which 
have been described here, clearly demonstrate the 
potential of Che technique for monitoring subtle 
earth movements. However, It Is equally clear 
that relevellng estimates of crustal movement are 
Influenced by near-surface movements and as yet 
poorly understood systematic errors which can 
obscure or be mistaken for tectonic deformation. 
Thus, uncritical interpretation of relevellng 
observations can lead to erroneous tectonic 
conclusions, which in the case of earthquake 
prediction could entail serious social 
ramifications. The checking techniques (e.g., 
circuit closure analysis) and reliability 
criteria illustrated in thi.3 study, represent an 
attempt to quantify specific procedures for 
evaluating the tectonic significance of 
particular leveling data sets. Although not 
foolproof, these procedures have proven effective 
in a number of cases at discriminating tectonic 
movements from suspect effects. However, even 
when spurious effects can be eliminated, relating 
observed deformation to preselsmlc mechanisms may 
be quite difficult because of our limited 
understanding of precursory phenomena and our 
general inability to distinguish them from 
vertical movements due to other causes (e.g., 
magmatic activity, isostatlc movements, aseismic 
erogenic deformation, etc.). Furthermore, the 
sparse distribution of leveling surveys in both 
space and time, even in areas like southern 
California, makes it highly unlikely that 
precursory movements for all but the largest 
earthquakes will ever be detected. In order for 
leveling to become more than an accidental 
contributor to earthquake prediction, a 
systematic leveling program designed for 
geodynamic rather than geodetic objectives is 
needed to develop the observational background 
required to recognize possible preselsmlc 
movement. 
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FIGURE CAPTIONS 


Figure !• Leveling routes for which crustal 
movement information has been obtained in the 
U*S« Locations of U«S* earthquakes for which 
releveling evidence of crustal movement has 
been reported are also shown. Numbers refer 
to Table 2. 

Figure 2. Magnitude of refraction error (Rg) 
normalized by height difference (AH) versus 
sight length (L) for various temperature 
differences. Based on relationship and 
theoretically determined constant given by 
Uoldahl (1980). 

Figuro 3. a) Profiles of apparent <*levation 
change and topography from Colorado Springs 
to Leadville, Colorado. Reversal of apparent 
tilt and correlation with topography strongly 
suggest elevation correlated error, 
b) Apparent elevation change versus elevation 
difference for 1954-1953 profile in Figure 
3a. Correlation coefficient (r) and 
regression slope (magnitude of error) are 
also shown. 

Figure 4. Elevation change profiles along 
east coast of U.S. (map at right). 
Unadjusted profile based on observed 
elevations from leveling assuming constant 
velocity iQovcmen:: (modified from Brovm, 

1978). Adjusted profile is same data 
adjusted by otlu-r levelings lying inland from 
coast (tide gauge data not used in 
adjustment). Squares show similar profile 
derived from tide gauges. 

Figure 6. Location of elevation change 
profiles which indicate subsidence possibly 
due to groundwater effects. Shaded areas 
represent previously published cases of near 
surface subsidence. 

Figure 7. Map of Los Angeles basin. Contours 
indicate depth to basement (M) , heavy dashed 
lines are faults, stippled areas are bedrock 
outcrops, heavy solid line (San Pedro to Los 
Angeles) is leveling route - crustal 
movements for period 1955-1964, Indicating 
sediment compaction, shown below map. 
Asterisk shows location of observation well 
for which waLuc level history is shown at 


right 


Figure 9. Contours of elevation change (I 
ran/yr and 5 inm/yr contours shown) for doming 
of Hebgcn Laku region* Movements may reflect 
preseismic movutaents Cor 1959 Hebgen Lake 
earthquake (kolllnger et al», 1977). 

Figure 10. Locations of leveling routes and 
benchmarks (dots) in Socorro-Albuquerque, New 
Mexico area. Outline of mid-crustal magma 
body is also shown (from Rinehart et al., 
1979). 

Figure 12. Leveling routes in Southern 
California for which crustal movements have 
been investigated in this study. 1971 San 
Fernando earthquake epicenter (*) and surface 
fault are shown along with contours of 
Palmdale Bulge as reported by Castle (1978). 

Figure 13. Relative movements for sequential 
time intervals and topography Cor route I in 
Figure 12. 

Figure 13. Relative movement for sequential 
time intervals and topography for route II in 
Figure 12. 

Figure 16. (Top) Leveling routes traversing 
the Saugus aquifer* 

(Bottom) Relative movements along routes 
shown at top. Note correlation between 
subsidences and aquifer location. Section 
labelled S.C.R. occurs within Che alluvial 
aquifer o£ the Santa Clara river and may be 
subject to subsidence due to compaction of 
this aquifer (Mark et al*. 1981). Section 
labelled RCF may reflect refraction error 
(Figure 14b and text). 

Figure 17. Plot of subsidence versus change 
in potentlomet r Ic surface (l.e., water level: 
AP) times aquifer thickness (T) for 
benchmarlcs in and immediately adjacent to 
aquifer (from Reillnger, 198C). Different 
symbols refer to different releveled 
segments: circles-Saugus to North 
(1953-1964); squares-Saugus to South 
(1955-1964); triangles-Saugus to East 
(1955-1961). The three circled points lie in 
the southeastern part of the aquifer and may 
reflect either different sediment 
characteristics or some effect ocher than 
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sediment coiapaccLon In this arei> 

Figure 18> Profiles of apparent elevation 
change from Llano to Azusa, California used 
to define the Palmdale Bulge (see Figure 12 
for location) • 

Figure 19* Correlation between apparent 
Qoveraent and terrain for segment (A to B) uf 
1962-1971 profile shown In Figure 18. 
Correlation coefficient (.98) and slope of 
regression line are also given. Correlation 
and reversal of apparent movement strongly 
suggest elevation-dependent systematic error. 

Figure 20. Relative movement for sequential 
time Intervals and topography for route III 
In Figure 12. 


Figure 5. Leveling loops used to investigate apparent tilting between Davis Junction, Illinois and 
Willard, Ohio* Elevation change profile and topography along Davis Junction to Willard route shown 
at right. Mlsclosures for each loop, given In brackets, suggest systematic error In 1967-1969 
leveling from Davis Jet* to Willard, Ohio. 
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Tublc 1. Magnitude of systematic discrepancies between repeated levelings (appar- 
ently due to leveling errors) in areas of subdued relief. Oiscrenancics are given 
in terms of apparent tilt (nun/km) and apparent tilt rate (10”®yr ^). Distance over 
which discrepancies accumulate monotonical ly and average relief along route are shown. 


AULA 

MAGNITUDI. 
(MM/ KM) 

MAGNlTUDi; 
(lO’O YU‘‘) 

U1ST/\NCE 

(K.M) 

AVLKACIi KId.lLE 
(M) 

East Coast 
U.S. 

O.S 

1.6 

650 

<5 

West Coast 
U.S. 

0.3 

3.3 

380 

<10 

Midcontinent 

U.S. 

1,0 

2.6 

i40 

<100 


Quoted Limit for Systematic Error ir a Single Survey <0.2 mm/km (Bomford. 1971) 


Table 2. Earthquakes for Which Vertical Movements Derived from Leveling have been Reported. 

Maximum Observed Deformation and Typical Dimension of Effected Area are also Given. Numbers refer to Figure 1 . 
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ELEVATION CHi\‘JCES NEAR THE SAN GABRIEL FAULT, SOUTHERN CALIFORNIA 

Robert Relllnger 

Department or Gejloglcal Sciences, Cornell University 
Ithaca, N.Y. U853 


Abstrac t . Analysis of repeated leveling ob- 
servations in the vicinity of the San Gabriel 
fault In Southern California indicate subsidence 
immediately south of the fault relative to points 
to the north, south, and east. These observa- 
tions were previously Interpreted as reflecting 
tectonic notions associated with either the 
"Palmdale Bulge" or with preselsmlc effects of 
the San Fernando earthquake. Relative subsidence 
between 1953 and 1964 reaches approximately 9 cm 
and extends over a distance of more than 20 km. 
Subsidence occurs directly above the Saugus aqui- 
fer and shows a temporal correlation with the 
history of water level decline within the aqui- 
fer. The degree of subsidence of individual 
benchmarks is roughly proportional to the product 
of aquifer thickness and water level decline at 
the location of the benchmarks. These observa- 


tions strongly suggest that movements of the 
surface near the San Gabriel Fault, previously 
inferred to be of tectonic origin, actually re- 
sult from near surface sediment compaction with- 
in the Saugus basin. 

Introduction 

Releveling estimates of crustal movement may 
reflect tectonic deformation, nontectonlc motion, 
or systematic measurement errors. Proper Inter- 
pretation of relcveling measurements entails de- 
scrlmlnatlng between these effects. The re- 
leveling measurements in Southern California 
which are presented here were initially inter- 
preted as representing tectonic movements asso- 
ciated with the "Palmdale Bulge" (Castle et al., 
1976). More recently. Strange (1980) has sug- 



Flg. 1. Contours (meters) showing thicknes-:; of the Saugus aquifer (kubson, 1972). Level- 
ing routes and benchmarks (dots) crossing tlje basin are also shown. Dashed lines are 
faults. 
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Fig. 2. A) Top: Movement of benthmarkB along 
approximately North-South route plotted ver»ua 
diutance along route. BMXl/0 on northern peri- 
phery of aquifer assumed stable. Years of level 
Ing surveya shovn at top. Middle" Profile of 
thlckneaa of Saugus aquifer along leveling -■ ate 
(from Robson, 1972). Bottom; Profile of decline 
In potentlometrlc surface (water level) of 
19b3 relative to 19^5 surface (from Robson, 
1972). B) Profiles along route from north of 
Saugus to Lang (see Fig. 1). Same format as A. 


gested that these measurements reflect prcselsmlc 
deformation associated with the 1971 S.sn Fernando 
earthquake. In this study, evidence Is presented 
which strongly suggests that movements of the 
ground surface south of the San Gabriel fault 
In Southern California result, for the most pa't, 
from near surface sediment compaction due to 
fluctuations of the water level within the Saugus 
aquifer and not from tectonic activity. Al- 
though these crustal movements most likely are 
not tectonic In origin, the do represent real 
surface movements and in this sense further de- 
monstrate the ability of precise leveling to 
detect relatively subtle deformation. 

Data 


(955 - ( 96 / 




Figure 1 shows contours of the thickness of 
the confined Saugus aquifer (Robson, 1972) with 
leveling lines and benchmarks superimposed. Pro- 
files of elevation change, aquifer thickness and 
estimated change In potentlometrlc surface 
(water level) along the survey routes are shown 
In Figure 2a, b. The elevation change profile 
from north of Castalc through Saugus Is plotted 
assuming stability for benchmark X3>0 Immediately 
adjacent to the aquifer. This benchmark Is 
chosen as a reference since It lies outside of 
the aquifer and appears generally stable rela- 
tt'/o to benchmarks further north. The elevation 
chinge profile from Saugus to Lang is plctted 
assuml'ig stabllit ■ near Lang since the bench- 
marks immediately adjacent to the Saugus aquifer 
occur within the alluvial aquifer of the Santa 
Clara River and appear to have subsided them- 
selves. The change In water level shown in 
Flgu re 2a, b represents the calculated decline as 
of 19'^3 relative to the steady state level es- 
timated to have existed around 1945 (Robson, 
197'). This decline was due to pumping from the 
Sa'i'ij? aquifer as well as overlying .'Pluvial 
aq'. - ors, and to an extended drought i. obson, 

19 72, i’p. 8, 40). Because the water level was 
alrt'iily below the steady state at the time of 
the first leveling survey (1953), the decline In 
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Fig. 3. Profile of elevation change fnr period 
146^ to !'*69 along route fi.'.ii nortli of Castalc 
to Saugus (see Fig. 1 for location). 
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-AP* T(cm*) . 10^ 

Fig. 4. Plot of subsidence versus change In 
potent locet rlc surface tines aquifer thickness 
for benchmarks In and Imaiedlately adjacent to 
aquifer. Different symbols refer Co different 
releveled segments: c Ire lea-north of t'astalc to 
Snugus (195 3-1964); square.s-Saugus tc South 
(1955-1964); t rlangles-Saugua Co bang (1955- 
1961). 


water level shown In Figure 2a, b Is sontwhat 
larger than the actual change between leveling 
surveys. However, water levels declined more or 
less continuously up to 1963 (Robson, 1972). 

Figure 2a, b Indicate a reasonable spatial 
correlation between subsidence, aquifer thickness 
and change In water level. In addition, the 
Viter level within the Saugus aquifer stabilized 
after 1963 and began to recover around 1968 
(Robson, 1972). Examination of post-1964 eleva- 
tion changes across the northern jiart of the 
b.Tsin Indicate a gener.nl absence of subsidence 
of the bsiln (Figure 3). Thus, the observed 
e’ vitlon changes correlate spati.slly with 
a," for geo-iecry and tr. porally w th the hl.story 
ot vnter lev*', decline. 

In C''e i ;e->l case, the amount of compression 
(.subs l.'.ence) of a sediment layer is Proportional 
t.) the thickness of the layer and the increase 
' f* . • : e pressure ( rain to grain ressure) 
o ’-ho s- -icnts (Terza.’.hl and Peck, 1967; p. 
?'-). : .erlng the waC' r level produces a 

; ; 'port >jdl rede 'ion In the buoyant force be- 
rw.'fa gralrn and .lencc cause: an Increase In the 
effective pressure on the u.: -■•■lying sediments 
(U-rraghl .i.'d Peck, 1967 ; p. 384). Assuming 
unlfo-.n ;e !*r,ent chrracti rlstlcs throughout the 
aq-.'er (-'cst.int coepr. islbllity) , to a first 
tlfcat i >n .he amount of compaction (subsl- 
d;:. -e) at 3 p.irtlcular location •..111 be roughly 
pr '.)0 r t ional to the product of water level de- 
c' -e .i".d aquifer thickness. Figure 4 gives a 
,’icr o' re', i-lz' subsidence versus the product 
of the cb < > in water level and nqcil.'’er thlck- 
r.-' •' '•onch.': „rks In and immedlac. ly adj.ncent 

to- ■ . •• 1 .S aquifer. Considering uncer- 

tJ... •' 1 1.. aquifer parameters and pr..b.tble 

sp-itlal variations In aquifer compressibility, 
the ^bs.‘r' ’.'t.ins fit the simple linear relatl.m- 
sh.p 'r-.- .'1. ■ • tt'.ree circled points on the 

r*..,>': s l-de ■ the ; ■ .it 1 .e In • •: .suutheas t e i •! 

1 ' ■ '.,.r :'i 1) and may reflect 

'.t.-ier different sediment cf ir.icterist Ics or 


possibly tectonic movement In this area. The 
slope of the straight line In Figure 4 Is an 
cstlm.'ltc ol the average compressibility of the 
Saugus aquifer. Although the compressibility 
can vary thiough at least two orders of magni- 
tude depending on sediment type, stress level, 
and stress history, the value obtained from 
Figure 4 (’7 x 10”'^cm"’) is quite comparable to 
those reported for other confined aquifers In 
California (Lofgren, 1979, p. 35). 


Cone luslons 


Repeated leveling surveys conducted between 
1953 and 1964 Indicate subsidence of the Saugus 
Basin relative to bcnchtmirks on Its periphery. 
Subsidence extends over a distance of about 20 km 
with maximum relative subsldeiue near the center 
of the basin reaching 9 cm. The spatial corre- 
lation between the area of relative subsidence 
and the Saugus aquifer, the temporal correlation 
with the history of water level decline, and the 
parameter dependence Indicated In Figure 4, 
strongly suggest th.it movements In the Saugus 
Basin, previously inferred to be a*jsociated with 
cither the "Palmdale Bulge" or presel,imlc effects 
of the 1971 San Fernando earthquake, result from 
near-.surf .ICC sediment compaction. This study 
demonstrates the need for caution when applying 
tectonic Interpretations to relevellng observa- 
tions In areas of unconsolidated sediments. On 
the other hand, the fact that the leveling meas- 
urements accurately detected real surface move- 
ments under actual field conditions demonstrates 
the potential of this technique for monitoring 
crustal deformation. 
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ABSTRACT 


Deformation possibly related to an earthquake in Texas has been 
revealed from analysis of precise leveling data collected by the National 
Geodetic Survey. Relative elevation cha >ges deduced from first-order 
leveling surveys conducted between 1910 and 1957 show that at least two 
areas in western Texas near the epicenter of the 1931 Valentine earthquake 
(M = 6.4) subsided relative to their surroundings. Apparent subsidence 
east of Van Horn, Texas is attributed primarily to the combined effects of 
groundwater withdrawal and topography-related survey errors. In contrast, 
relative subsidence near Valentine, Texas, of 11.2 + 1.0 cm extending over 
a distance of about 20 km does not appear to be due to either near-surface 
effects or leveling errors and thus may represent coselsmic deformation of 
the Valentine earthquake. The predominance of subsidence over uplift 
indicated by the leveling data suggests a normal fault mechanism for the 
earthquake. This Is consistent with the Cenozolc structure of the region 
which is characterized by Basin and Range morphology. If the observed 
subsidence was caused by the 1931 earthquake, the leveling data suggest 
that the epicenter for this event lies considerably closer to Valentine 
than originally reported. This new location is supported by a recent 
relocation for the Valentine earthquake, and by local intensity reports. 
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INTRODUCTION 

Western Texas, (Figure 1) structurally a part of the Basin and Range 
Province, is considered by some to be the eastern offset of the Rio Grande 
Rift (Woodward et al., 1975; Muehlberger et al., 1978). Faults of Tertiary 
and Quaternary age shape the landscape into alternating horsts and grabens. 
The grabens are commonly filled with large, thicknesses of unconsolidated 
alluvial, lacustrine clay, silt. sand and gravel. These deposits are the 
major sources of groundwater in western Texas. Eviience of recent tectonic 
activity in this area includes north- and nor thwest-trend ing faults 
offsetting Quaternary georaorphic surfaces (King, 1965; Muehlberger et al., 
1978), historic seismicity (Coffman and von Hake, 1973; Sanford and 
Toppozada, 197A), and recent uplift in the Diablo Plateau area (Brown et 
al., 1978; Reillnger et al., 1980). 

Since very little is known about the long-term earthquake history of 
v'estern Texas, it is difficult to evaluate seismic risk from analysis of 
earthquake statistics. The Valentine, Texas, earthquake of August 16, 1931 
(M = 6.4), was the strongest reported earthquake to occur in the region 

(Sellards, 1932). Byerly (1934) located the epicenter at 30.9*N and 

o 

104.2 W, in the Davis Mountains, 38 km northeast of Valentine, Texas 
(Figure 1). More recently, Dunas et al. (1980) have used regional ard 

local teleseismic P-wave arrival time anomalies from the Gnome underground 
nuclear explosion to determine station corrections and relocate the 1931 
event at 30.69 N, 104.57 W (Figure 1). Tliis relocation is considerably 
closer to Valentine than the previously determined epicenter. First P 
motions analyzed by Byerly (1934) and later reexamined by Sanford and 

o 

Toppozada (1974) suggest predominantly normal faulting, striking N40 W and 
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o 

dipping 74 SW. Dumas et al . (1980), using the same first P motion 

observations, have determined a new fault plane solution and suggest that 

o 

the earthquake occurred on a predominantly strike-slip fault striking N'i9 W 

o 

and dipping 70 NE. 

The incon'i isteocy of faul. planes determined by seismic data and the 
lack of evidence of surface faulting in the Valentine region make it 
difficult to assign this event to a specific fault. A very poorly located 
earthquake on 1 August 1975 (M=»' .8), was felt in Valentine, however m 
focal mechanism was determined for this event. Nevertheless, the Valentine 
region is seismically active at present (Dumas et al . , 1980). 

In this study we present elevation changes deduced from precise 
leveling surveys in the vicinity of the Valentine earthquake. Because 

leveling measurements can be influenced by near surface movement ^ 
systematic errors, and tectonic deformation, specific criteria are applied 
to separate these effects. Anomalous relative subsidence near Valentine, 
Texas suggests the existence of an active burled fault that may have been 
the source of the 1931 Valentine earthquake. If so, this is the first 
documented case of contemporary earthquake deformation in Texas. 


DATA ANALYSIS 

The primary data used in this study are estimates of crustal movement 
derived from repeated leveling surveys conducted by the National Geodetic 
Survey and water-well levels reported by the U.S. Geological Survey and the 
Texas Water Comniosion. 

The details of the leveling procedure and the method used to derive 
crustal movement information from repeated leveling surveys have been 
described elsewhere (Bomford, 1971; Broi/ii and Oliver, 19'^6). Briefly, 
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leveling Is an optical technique used to determine the differences in 
elevation between closely spaced ( 1 km) points (benchmarks) on the 
earth's surface. Repeating a survey at some later time (^t =. tens of 
years) results in new elevation observations that can be compared to the 
original observations to estimate relative movement. 

Figure 1 shows the locations of the leveling routes used In this 
study, and Figures 2a and 3a', elevation changes along these routes. All of 
the survc/s used to determine elevation changes were performed between 1910 
and 1957 by the National Geodetic Survey according to first-order stan- 
dards. There is no evidence of abnormal benchmark instability in their 
published descriptions or from field inspection ( J . Dorman, 1979, Personal 
communication) . The observed elevation changes in Figures 2a and 3a have 
been plotted with respect to benchmark Y17 (Sierra Blanca) as a reference. 
This does not restrict the data since only relative movements are 
s Ignif leant . 

Apparent crustal movements between Sierra Blanca and Sanderson, Texas, 
are shown by the profile in Figure 2a, which reveals, among other things, a 
broad regional down-tilting to the southeast. The magnitude of this tilt 
is about 0.76yurad. It is difficult to ascertain whether this tilting is 
due to systematic leveling errors or to a regional tectonic effect (Brown 
et al . , 1980). For example, since the route is a partially north-south 
line, tidal and unequal lighting errors (Bomford, 1971; Balazs, 1975) could 
have some influence on the data. In either case, this signal is almost 
certainly not directly related to the Valentine earthquake. Because of 
this, the data were adjusted to remove the regional tilt. Figure 2a shows 
the adjusted elevation change profile along with the unadjusted profile. 
Both profiles indicate anomalous subsidence near Valentine, Texas, although 


PAGE 6 


with slightly differont magnitudes. The significance of this feature and 
its implication *or the 1931 Valentine earthquake will be discussed in the 
following section. 

The profile of movements between Sierra Blanca and Monahans (Figure 
Ja) shows no regional tilting, possibly a consequence of its predominantly 
east-west orientation. However, this profile does show a negative 
correlation with topography between benchmark BM H17 (near Van Horn) and BM 
S16 (near Pecos). There are a number of known systematic leveling errors 
that accumulate with elevation (unequal atmospheric refraction; rod 
miscalibration ; Boraford, 1971). Although a correlation with topography 
does not in itself invalidate tVie leveling observations, it does constitute 
a negative factor as far as credibility is concerned (Brown et al., 1980), 
In an attempt to remove the effects of possible systematic errors, a plot 
of relative movement versus elevation was constructed for the segment of 
the movement profile that correlates with topography (Figure 4). The slope 
of the regression line fit to this plot is 25 mra/100 m, corresponding to 
the average magnitude of the elevation-correlation. Removing the part of 
the signal wl.lch correlates with topography results in the adjusted profile 
segment shown in Figure 3^. The broad subsidence features between BM H17 
and BM S16 is no longer apparent. The adjusted profile shows few 
outstanding features except for a small relative subsidence between BM 3826 
and BM H17, and a large relative subsidence centered near Pecos, Texas. 
Subsidence near Pecos has been shown to be the effect of sediment 
compaction due to groundwater withdrawal (Rosepilcr and Reilinger, 1977). 
Subsidence between BM 3826 and BM H17 also appears to be due to near 
surface movements and will be discussed further in the next section. 

The other data used in this study are water levels in wells located 
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near the leveling routoa. The geographical distribution and thickness of 
the principal aquifers In the vicinity of Van Horn and Lobo, Texas, are 
given by Gates et al (1978) and Hood and Scalaplno (1951) . Pumping for 
Irrigation from these aquifers began in 1949 (Gates et al., 1978), 
Increased rapidly In the early 1950's and has remained fairly stable to the 
present. As water withdrawal Is larger than natural recharge, well levels 
have declined in these basins. Water level declines reached 15 meters 
between 1949 and 1957 near Lobo, Texas (Figure 2b). The decline In water 
levels for wells located near the leveling routes are shown in Figures 2b 
and 3b (data from Gates et al. (1978) and Hood and Scalaplno (1951)). 

DISCUSSION AND INTERPRETATION 

Elevation C hange Due to W ater Level Dec line 

Vihencver the water table Is lowered, the effective pressure on the 
u: derlying sediments is increased. The increase in effective pressure, or 
grain-to-grain loading, causes the sediments to compact which results in 
subsidence of the land surface. 

Land subsidence due to withdrawal of groundwater has been identified 
in the Pecos region (Rosepiler and Rellinger, 1977) and in various other 
parts of North America where large volumes of water are being removed from 
underground aquifers (Poland and Davis. 1959). Therefore, care must be 
taken to avoid misinterpreting subsidence due to fluid withdrawal as 
tectonic notion. Inspection of Figure 2b shows a strong correlatlo i 
between water-level decline and land subsidence between BM L23 and BM Q23. 
This strongly suggests that the land subsidence in the Lobo Flat region is 
due to compaction of the aquifer in response to decline in hydrostatic 
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pressure. Compaction probably occurs chiefly in the intrabeddcd clay and 
reworked volcanic elastics. The ratio of land subsidence to decline in 
water level is about 1 to 400, equivalent to 25 mm subsidence for 10 m of 
water-level decline, a value comparable to that observed in other areas in 
the United. States (Poland and Davis, 1969). It is important to note that 
the apparent subsidence south of Lobo Flat, i.e., near Valentine, is not 
correlated with water-level decline and therefore is most likely not due to 
pumping (Figure 2b). 

A less-extensive land subsidence due to water withdrawal is indicated 
west of Van Horn along the Sierra Blanca-Pecos profile (Figure 3b) by a 
similar correlation of subsidence and water-level decline. The primary 
effect is confined between BM 3826 and BM H17. The ratio of land 
subsidence to water-level decline at BM 3826 Is 1:700. 

Implications for the 1931 Valentine Earthquake 

As noted in the previous section, the anomalous relative subsidence 
centered near Valentine is not easily attributed to near-surface effects, 
and the magnitude of this subsidence is considerably greater than random 
error. It is therefore suggested that the observed relative subsidence 
represents tectonic deformation associated with the Valentine (;arthquake. 

The predominance of subsidence over uplift indicated by the leveling data 
suggests a normal fault mechanism for the earthquake (Savage and Hastie ,1966) . 

Both fault plane solutions for the 1931 Valentine earthquake (Sanford 
and Topozada, 1974; Dumas et al., 1980) suggest a Northwest striking fault. 
This trend is in good agreement with both the regional structural grain and 
the long axis of the isoseismal contours for the 1931 earthquake given by 
Sellards (1932). The leveling route traversing the Valentine area is also 


oriented in a 
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generally northwest direction. Therefore, If the observed subsidence is 
due to fault slip as sugi^ested, the areal extent of subsidence (^30 km) 
is a rough estimate of the length of the earthquake fault. Kanamori and 
Anderson (1975) give an empirical relation between fault dimi>>' >lons and 
surface wave magnitudes. Using these relationships and assur a fault 
width of 10 km suggests a magnitude in the range 6.3 - 9 for the 
Valentine earthquake. Although not unreasonable, this is somewhat larger 
than the magnitude estimated for this event (Mb^S.S to 6.4, Cutenburg and 
Richter, 1949; Dumas et al., 1980). However, the fault length determined 
from the zone of subsidence is probably overestimated for the following 
reasons: 1) surface deformation will extend beyond the ends of the fault; 
and 2) observed subsidence could include the effects of preselsralc and/or 
postseismic slip on the fault or on extensions of the fault. 

Although the details of the faulting responsible for the Valentine 
earthquake arc poorly constrained by the available data, the observed 
pattern of surface deformation is at least consistent with a northwest- 
striking normal fault dipping southwest and located just north of 
Valentine, Texas. This location is in agreement with the area of maximum 
intensity for this event (Sellards, 1932), and with the recomputed 
epicenter reported by Dumas et al. (1980). 

CONCLUSIONS 

Relative elevation changes deduced from repeated leveling surveys in 
the Van Horn-Valentine region of western Trans-Pecos, Texas, indicate 
significant vertical movements between 1910 and 1957. Although some of the 
movements appear to be due to water withdrawal (e.g.. 


rear Van Horn and 
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Lobo, Tftxaa) , the most anomalous movements occur near Valentine and do not 
appear to be associated with groundwater withdrawal or leveling errors. 
Subsidence near Valentine, reaching 11.2 + 1.0 cm relative to points to the 
•'orthwest and southeast. Is thus most likely of tectonic origin. The 
pattern of deformation near Valentine Is consistent with normal faulting on 
a northwest-striking, southwest-dipping fault. In agreement with one 
interpretation of the focal mechanism inferred for this event. This 
earthquake Is the largest recorded in Texas and the observations presented 
here may represent the first 'ase of earthquake deformation in this region. 
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FIGURE CAPTIONS 


Figure 1. Regional geography and Cenozolc faults of the study area. 
Heavy lines Indicate leveling routes. Open squares indicate cities. 
Crosses represent location of benchmarks. Epicenter locations of the 1931 
Valentine, Texas, earthquake reported by Byerly (I97^ » J'nd Dumas ot al. 
are also .‘epresented. 

Figure 2. (a) Profiles of relative vertical crustal movements from 

Sierra Blanca to Sanderson, Texas, derived from N.G.S. surveys conducted in 
1917 and 1957. Unadjusted profile is shown by crosses; adjusted profile, 
by solid circles. The upper curve reflects the topography a’ong the 
profile. Locations of some cities and benchmarks discussed in the text are 
indicated, (b) Expanded display of recent crustal movements and water- 
level decline (open circles) in the vicinity of Lobo , Texas. 

Figure 3. (a) Profiles of relative vertical crustal movements between 

Sierra Blanca and Monacans, Texas, derived from N.G.S. surveys conducted in 
1910 and 1956. Topography is shown between the unadjusted profile 
(crosses) and adjusted profile (solid circles). Location of some cities 
and benchmarks discussed in the text are indicated. (b) Expanded profile 
of recent crustal movements and water-level decline (open circles) near Van 
Horn, Texas. 


Figure 4. Plot of net elevation change versus elevation. 
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DETECTION OF VERTICAL FAULT DISPLACEMENTS BY PRECISE LEVELING: 
A CASE STUDY IN WESTERN KENTUCKY* 


ABSTRACT 


Relative vertical displacements of bench marks in extreme western 
Kentucky have beer determined by comparison of successive leveling 
surveys in 1947 am 1968* The resulting pattern of apparent surface 
deformation shows .< steep offset which can be closely modeled by a 
nomal fault burled in an elastic half-space. The offset is located 
near the northern boundary of the Misslsslpoi Embayment and the New 
Madrid seismic zone, an area where faults have previously been inferred 
on the basis of both geological and geophysical evidence. If the 
apparent movement is due to slip along a fault, several lines of 
evidence (regional structure, earthquake data, and llneations) suggest 
that the postulated fault trends north-northeast. Tliirteen earthquakes 
were recorded in this area between the times of leveling; focal 
mechanisms exist for three of these. The nearest of these three focal 
mechanisms to the leveling offset implies normal faulting. The 
magnitude of the earthquake, however, appears to be too small to 
account for the amount of slip required by the fault model. Thus the 
apparent deformation may have accumulated with several undetected small 


* By F. Steve Schllt and Robert E. Rellinger 
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earthquakoB, or gradually as aBslsmic creep 


INTRODUCTION 

The purpose of this study is to examine the possibility that 
anomalous surface deformation near the northern end of the Mississippi 
Embayment, as measured by repeated precise leveling surveys » was caused 
by movement along a burled fault* The area studied is between 
Wickliffe and Paducah, Kentucky, and the leveling data, when 
supplemented with other geological and geophysical information, suggest 
dip-slip displacement along an unmapped, concealed fault. 

Brown and Oliver (1976) have compared regional structural features 
of the eastern United States with relative rates of vertical crustal 
movement as determined from multiple first-order leveling surveys. 
They found a significant correlation between the velocity patterns and 
large-scale (hundreds of kilometers) geological structure. A portion 
of their data for the central United States is examined here on a more 
local scale, with Interest in shorter wavelength behavior. 

Ccselsmic surface deformation has been studied using leveling and 
other techniques for many large, dip-slip earthquakes (e.g.. Savage 
and Uastle, 1966; Rcilinger and Brown, 1981). In most of these 
examples, movements are large and are easily Identified by their close 
temporal and spatial association with earthquakes, and in some cases, 
surface faulting. This study Involves a more subtle degree of 
deformation and many kinds of evidence are considered in order to 
determine if the deformation is tectonic and if so, how it occurred. 


T 
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REGIONAL STRUCTURE AND GEOLOGY 

The MlsBlBBippi Embaymcnt (Figure 2.1) began to fora in the Late 
Cretaceous WIten the Fascola Arch (then connecting the Ozark and 
Nashville domes) was downwarped. The cause of this subsidence is not 
well understood, but it has been viewed as an Isostatic adjustment 
subsequent to Precaiubrlan rifting, which extended into the continent as 
the failed arm of a rrij. -e-junction (Burke and Dewey, 1973; Ervin and 
McGinnis, 1975). Depression continued as a linear trough until late 
Eocene, with an accumulation of about 1000 meters of sediment in the 
central portion of the embayment. A cover of alluvial deposits is Late 
Tertiary in age, and is locally overlain by Pleistocene loess. 

Nearly all of the observed surface faults in this region occur in 
the Paleozoic rocks surrounding the Cretaceous and younger fill of the 
Mississippi Embayment (Ueyl, 1965; Hook, 1974). The major 
displacements on these Paleozoic faults probably occurred in 
Pennsylvanian and Mlssisslppian time (King, 1969). Far fewer faults 
are mapped within the Cretaceous and younger rocks of the embayment 
(York ar.fi Oliver, 1976; Russ, 1979; Zoback, 1979). Recent seismicity 
(Figure 2.1) defines several linear NE-trending zones approximately 
parallel to the embayment axis, and some zones trending NW (Stauder and 
others, 1976; Herrmann and Canas, 1978). 

The difficulty of finding mappable faults which can be clearly 
related to this seismicity may be largely due to a masking effect of 
the unconsolidated alluvial cover. In an attempt to map faulting 
within the embci^ment, Fisk (1944) made a detailed study of llneations, 
and obtained a grid-llkc pattern of orthogonal faults for the entire 
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figure 2.1. Selsmotec tonic map of the northern Mississippi 
embayment and surrounding regions (Hlldenbrand and others, 1978). Study 
area is indicated by rectangle. 






lowur MiswiBBippi River Valley. Tlie overall pattern of Fisk (1944) in 
not ref lec Ceil In the eeiemlclcy* and more recently Stearns and Zurawski 
(1975) have suggested that lineatlons may have often been Interpreted 
us faults when they are actually Joints or Joint swarms. Stearns and 
Zuratmki (1975) therefore used subsurface structural contours based on 
well data in conjunction with lineations, and their inferred fault 
pattern is considerably different from that of Fisk (1944), and shows 
somewhat better correlation with seismicity. 

York and Oliver (1976) reviewed the evidence for Cretaceous and 
Cenezolc faulting in the central United States, and found all known 
faults of this age to lie in the northern tip of the eobayroent. The 
inferred fault pattern of Stearns and Zurawski (1975) likewise 
indicates numerous faults Just inside the northern margin of the 
embayiucnt. Zuba'.k (1979) delineated two fault systems by seismic 
reflection profiling near Rcelfoot Lake in northwestern Tennessee, 
about 75 km southeast ct our study. Thus, the leveling data in this 
region deserve careful examination for evidence of modern fault 
activity. 


LOCAL STRUCTURE 


The dense Paleozoic faulting of the western Kentucky faulted area 
disappears at the northern edge of the embayment (Figures 2.1 and 2.2). 
This area has been mapped in detail and shown to be a northeast 
trending horst and graben fault complex with vertical displacements of 
150-500 meters, associated with northwest-southeast lateral extension 
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o£ opproximately l«5 km (Uook« 197A)» Normal i’’aultin)} (with dip» 
averaging 70-75 degrees) ia predominant, although minor atrike-allp 
motion ia evident* The fault aone is now apparently quiescent, aa 
almost no seismic activity has been detected there* 

Considerable evidence, however, suggests the continuation of 
basement faulting beneath ond/or into the embayment* In southernmost 
Illinois, subsurface mapping has revealed a complex of buried grabens 
and faults just within the northern boundary of the embayment* At 
least a past of these displacements occurred after the formation of the 
sub-Cretaceous erosional surface, with episodic activity in the Late 
Cretaceous, post Eocene, Pliocene, and Pleistocene (Ross, 1963a; Ross, 
1963b) * 


SEISMOTECTONICS 


The pattern of seismicity from 1974 to 1976 is shown in Figure 2*1 
(Ulldenbrand and others, 1978; Stauder and others, 1976)* The major 
characteristics are the broadly-defined northeast trend within the 
embayment and a more diffuse pattern outside it* The existence of 
emaller-scale northeast- and northwest-trending lineaments is evident* 
Street and others (197<,) indicate hypo cen. cal depths for earthquakes in 
the New Madrid seismic zone ranging from 5 to 38 km, but the majority 
of the events shown in Figure 2*1 are probably between 2 and 10 km deep 
(Stauder and others, 1976). 

The inset rectangle seen in ' Figure 2*1 indicates the locality for 
which leveling data will be presented (Figure 2*2)* The center of the 
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titudy ucf.u ib approxtuatcly 75 ku northeast of New Madrid, Missouri, 
the situ of the large and damaging earthquakes of 1811-1812. 

Figure 2.2 shows the seismicity in the study area for two 

different time periods. Solid circles represent events which occurred 

between 1947.5 and 1968.7 (listed in Table 2.1a), the dates of first- 

order leveling surveys along the route between Wickliffe and Paducah, 

Kentucky. Open circles represent events for 1928*1947.5 and 1968.7- 

1972 (listed In Table 2.1b), time periods outside the leveling 

interval. If non-instrumental historical seismicity were Included in 

Figure 2.2, the greatest activity would be at Cairo, Illinois (37.0*' N 
o 

lat, 89.8 W long). For 1855-1972, ?9 events have been reported at 

Cairo, including one d-tmaging event in 1883 (intensity VI -VIII). The 

sharp clustering at Cairo since 1855, however, may reflect population 
bias, in absence of instrumental locations before 1928. 

Although thrust, normal, and strike-slip focal mechanisms are all 
present In the area shown by Figure 2.2, six of the eight nodal planes 
(four focal mechanisms) strike north-northeast between 356 and 031, a 
span of 35 degrees. The majority of nearby focal mechanisms (within 
about 100 km of Figure 2.2) Indicate thrusting, as noted by Street and 
others (1974). If the strike-slip event in southern Illinois can be 
associated with movement along tie trend of the New Madrid fault zone, 
the sense of diplaccment is consistent with east-west compression. 
Composite focal plane solutions in Arkansas and Missouri show a 
significant compomeut of right-lateral fault motion indicative of EW 
compression in the central part of the embayment (Herman and Canas, 
1978). The stress pattern indicated is not simple, but there is 
apparently a significant component of ESE-WNW compression in the 
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TABLE 2* la* Earthquakes in area near Cairo, Illinois (Figure 3) 
Time interval: 1947-1968. 7 


Date 

(yr .mo .da) 

Lat 

(N) 

Lon 

(W) 

Depth 

(km) 

Inten- 

sity 

Magni- 

tude 

Refer- 

ence 

1947.01.16 

36 59 

89 11 


II-III 

(3.2) 

a,b 

1953.05.06 

36 59 

89 11 

- 

III 

(3.4) 

a,b 

1953.05.15 

36 59 

89 11 

- 

III 

(3.4) 

a,b 

1957.03.26 

37 05 

88 36 

- 

IV 

(3.8) 

a ,b 

1958.01.27 

37 03 

89 12 

- 

V 

(4.2) 

a,b 

1962. U2. 16 

37.0 

88.7 

25 

I'l-IV 

(3.6) 

b.d 

J963.01.31 

36.9 

89.0 

- 

- 

3.0 

b,c 

1963.05.02 

36.7 

89.4 

- 

- 

3.1 

a ,b 

1963.08.02 (F2) 

37.0 

86.7 

18 

V 

4.0 

b,c,d 

1965.08.13 

37 19 

89 28 

33 

IV 

3.2 

41 ,b ,<l 

1965.08.14 (F3) 

37.2 

69.3 

38 

VTI 

3.8 

b.c.d 

1965.08.15 

37 22 

89 28 

16 

V 

3.5 

a,b ,d 

1965.08.15 

37 24 

89 28 

*• 

V 

3.4 

a ,b 


(FI): focal mechanism of Street and others (1974) (see Figure 3) 

(3.2): magnitude estlmarsd by Nuttll (1974) 

a: Docekal (1970) 

b: Nuttll (1974) 

c: Street and others (1974) 

di Stearns and Wilson (1972) 


TABLE 2.1b. Earthquakes in area near Cairo, Illinois (Figure 3) 
Time interval; 1928-1947 and 1968.7-1972 


Date 

(yr.mo.da) 

Lat 

(N) 

Ix>n 

(W) 

Depth 

(Ion) 

Inten- 

sity 

Magni- 

tude 

Refer- 

ence 

1928.04.15 

37 17 

89 32 


IV 

(3.8) 

a,b 

1930.08.29 

37.9 

89.0 


V 

» 

d 

1930.09.03 

36 58 

88 54 


III 

(3.4) 

a,b 

1930.09.03 

36 58 

88 54 


III 

(3.4) 

a,b 

1931.04.06 

36 51 

89 01 


IV 

(3.8) 

a,b 

1933.10.24 

37 17 

89 32 


III 

(3.4) 

a,b 

1934.08.19 

36 56 

89 12 


VI 

(4.7) 

a,b 

1934.08.19 

36 59 

89 09 


II-III 

(3.2) 

a,b 

1936.08.02 

36.7 

89.0 


III 

(4.1) 

a ,b 

1936.12.20 

37 17 

89 32 


11 

(3.0) 

a,b 

1939.04.15 

36 48 

89 24 


ITT 

(3.4) 

a,b 

1940.02.04 

37 13 

89 29 


III 

(3.4) 

a,b 

1940.05.31 

37 05 

88 36 


V 

(4.2) 

a b 

1940.10.10 

36.8 

89.2 


II-III 

(3.2) 

a,b 

1941.10.21 

36 58 

89 08 


IV 

;3.8) 

a,b 

1941.11.22 

37 17 

89 32 


II-III 

(3.2) 

a,b 

1942.08.31 

36 59 

89 11 


IV 

(3.8) 

a,b 

1942.11.13 

36 59 

89 11 


IV 

(3.W— 

a,b 

1970.12.24 

36.7 

89.5 


IV 

3.6 

b 

19r:.06.18 (F4) 

37.0 

89.1 


III 

3.2 

c 


see TABLE lA. for explanation 


cmbaytuenl: propc-. in concrast to the Paleozoic SE-WW extenbion seen In 
the adjacent weutevn Kentucky faulted area* 


LEVELING DATA 

The route between Paducah and Wickllffe (Figure 2.2) has been 
surveyed twice by the National Geodetic Survey using first order 
leveling. The first survey was made in 1947 and the second in 196C. 
Elevation changes computed from these leveling surveys (relative to 
Paducah) are shown in Figure 2.3. Also shown is the topography along 
the route. The standard error for the apparent elevation change (Brown 
and Oliver, 1976) at Wickliff relative to Paducah is + 15 mm . About 
60 percent of the bench marks along this route are in concrete posts. 

The most striking feature in Figure 2.3 is the steep gradient 
between benchmarks A and B (31 + 6 mm in 9 km) corresponu. to a net 
tilt of 3.4 + 0.7 X 10-6 rad. Ihis tilt is about an order of magnitude 
larger than most regional tilts in the eastern United States (Brown and 
Oliver, i976), and is definitely larger than random measurement error. 


DISCUSSION AND INTERPRETATION 

The segment of the route between A and B lies between the trends 
of the Big Creek fault zone (Fisk, 1944) and the western Kentucky 
faulted area (Hook, 1974) (Figure 2.2). It is possible, in light of the 
evidence presented thus far, that the movement indicated by Figure 2.3 


RELATIVE MOVEMENT, MM 



ELEVATION, M 






iB faulL rulattid* However, possible errors and various alternative 
explanations for the leveling observations must be considered. 

Since both leveling surveys were double run (elevation differences 
between successive benchmarks arc measured twice for each survey), it 
is highly unlikely that the observed tilting is due to blunders 
(mistakes by the field party). The tilting occurs over too short a 
distance and in too flr.t .u. area to be attributed to many of the 
possibile systematic sources A error (c.g., unequal refraction, tidal 
effects), especially considering the tendency for such errors to cancel 
each other when computing the elevation change. 

Near-surface non-tectonlc processes also appear Inadequate to 
explain the observed tilting. There is no evidence of significant 
fluid withdrawal large enough to cause sedimenc compaction in the area. 
Oil or gas wells are not present in this part of Kentucky, and the 
observed elevutiu;' changes show no relationship to cither water-level 
fluctuation (seasonal or long-term) or aquifer geometry (Davis and 
others, 1973). The region studied is not karstland, arguing against 
surface slumping due to cavern collapse (Gerlach, 1970). There is 
little correlation between relative movements and type of monumentation 
used along this line, arguing against significant syst<.matlc trends due 
to benchmark instability. Crustal loading (e.g. reservoirs) has been 
shown to cause appreciable downwarps, but no such cases of large 
artificial loads were found near the offset. 

If non-tectonic mechanisms are inadequate to explain the observed 
elevation changes, as seems to be the case, it is reasonable to 
hypothesize that they may be due to fault movement. The sense of 
surface deformation (down to the west) would be consistent with: 


(1) a wcBtward dipping; notmal fault 

(2) an eastward dipping thrust fault, or 

(3) a stike-slip fault with a significant component of dip-slip 
motion of type U) or (2) 

The closest for:al mechanism (F4 in Figure 2.2; Street and others, 1974) 
is about 10 km southwest of bench mark A. It is a thrust type , but 
that earthquake occurred four years after the second leveling survey. 
Tlie next nearest mechanism (FI in Figure 2.2; Street and others, 1974) 
is about 20 km south of benctunark A. It is a normal type, and occurred 
within the time interval between surveys. It is therefore difficult to 
associate the offset seen In the leveling data with one or the other 
type of fault. In each case, however, one of the nodal planes would be 
consistent with the sense of inferred surface deformation. 

From a single leveling route it is not possible to determine the 
nrcnl treriJ of surface 'leform.itiun. Iluwuvur, the trends of the iiudnl 
planes of the focal mechanisms, the trend of the structural trough 
formed by the top of the Paleozoic, Che trend of the Mississippi River 
as it bounds wo'^itern Kentucky, and the trend of lineations defined by 
local stream patterns (Figure 2.4) are all UNE. Titus it seems likely 
that if the observed movements are due to fault activity, the 
responsible fault strikes NNE. 


69*00’ W 


68* 56' 


FIGURE 2.4. NNE-trendlng lineaments defined by stream pattern near 
bench marks *'A" aiid ”B" (other \>ench marks are shown by "X") . Streams 
were traced froatU. S. Geological Survey 7 1/2' topographic quadrangle 
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THEORETICAL FAULT MODEL 

Savage and llaatlc (I960) have given a convenient formulation for 
computing the surface deformation due to a dip-elip dislocation in an 
elastic half-space* By experimenting with various fault parameters, 
the pattern of apparent deformation was modelled with moderate slip 
along a westwar'i dipping normal fault. Before making comparasons of 
the observed and modelled displacements, the bench mark positions were 
projected onto the profile P-P' (Figures 2.2 and 2.5), on th^ 
assumption that the postulated fault trends at 015. 

Tlie pattern of deformation was best matched by normal faulting, 
although thrust faulting can also give a plausible fit. Two normal 
fault models which fit the data reasonably well are described and shown 
in Figure 2.5. In addition to the problem of non-unlquencss, there is 
the ambiguity introduced by the lack of data extending further to the 
west. However, with certain assumptions, order of magnitude estimates 
of fault parameters are possible. Assuming that the fault is in the 
upper crust and that its length is a few times its width, a slip of the 
order of 10 cm over an area of 100 km2 is required to match the 
magnitude of the observed movements. 

If the deformation represents the strain release for one 
earthquake, the seismic moments implied by the two models in Figure 2.5 
are 6 X 10+24 dyne-cm and 1 X 10+25 dyne-cm. Applying an empirical 
relationship between seismic moment and magnitude, a magnitude of 5.5- 
6.5 would be indicated (Kanamorl and Anderson, 1975). Tl\e largest 
reported earthquake for this region between the times of leveling had 
an estimated magnitude of only 4.2. Thus, if the observations reflect 
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FIGURE 2.5. Suriacf? displacemenc predicted by dip-slip on a buried 
normal fault. Dots Indicate relative movements of bench marks 
(projected along the Inferred fault strike shovra in Figure 2.2) as 
measured by the leveling surveys . Displacement curves for two different 
normal fault models are shown. These two models give some idea of the 
tradeoffs possible in obtaining a reasonable fit to the observed 
movements . 
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fault notion, the dcfotiuatian must have accumulated niulli|>lu 

Buull events or as aseiautic cretta* 


CONCLUSIONS 

Elevation changes estimated from repeated leveling in western 
Kentucky may represent subsurface faulting. This conclusion, however, 
is based on a relatively limited data set, and the possibility of near- 
surface effects cannot be totally ruled out. 

Modeling of a burled fault suggests a dip-slip of about 10 cm over 

X 

a fault area of about 100 km . A normal fault is preferred, but 
thrusting is also possible. Tlie equlval(mt magnitude for a single 
earthquake is larger than any earthquake experienced during the 
Interval between surveys, suggest5.ng Incremental or aselsmlc 
displacement. It is recommended that these benchmarks be resurveyed as 
soon as possible to sec if further changes have occurred since 1968. 
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